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Abstract 

00 ' We obtain asymptotic expressions for the Green kernels of certain non-translation invariant 

, transition matrices using methods of semiclassical and microlocal analysis. Combined with 

a result by Bach and M0ller this yields asymptotic formulas for the truncated two-point 
correlation functions of certain non-translation invariant lattice models of real-valued spins. 
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1 Introduction 

' This article is a further contribution to the study of correlation asymptotics for lattice systems of 

, real- valued spins via the Witten-Laplacian approach pl[3ll4llT3l[20l[2T]. Like the authors did in the 

\^ ' cited articles, we consider a lattice system at very small temperatures and assume that the on-site 

, potentials have only one non-degenerated minimum and that the ferromagnetic pair interaction 

On ' is sufficiently weak. In this situation one expects the truncated two-point correlation functions of 

, the associated Gibbs measure to decay exponentially. In fact, in [2 [20] the leading asymptotics, as 

' the temperature tends to zero and the distance on the lattice tends to infinity, of the exponential 

decay were computed under similar hypotheses. In |21j even a full asymptotic expansion has been 
obtained. The results of [H [501 HI] hold for translation invariant systems on finite discrete tori and 
are uniform in the number of spins. (In [5T] the thermodynamic limit of the derived expansions 
' is discussed, too.) They provide a precise description of the Ornstein-Zernike behaviour which is 

, actually well-known for translation invariant systems since a long time; see, e.g., |17j . The starting 

point of the analysis of the correlation in [21 |3l H] [201 EI] is the Helffer-Sjostrand formula [71 [20] 
which represents the truncated two-point correlation functions as matrix elements of the resolvent 
of a certain Witten-Laplacian on one- forms. If one works directly in the thermodynamic limit, 
or more precisely, if one considers a spin system on the lattice T'^ whose equilibrium distribution 
^ ' is given by a translation invariant tempered Gibbs measure, then the formula for the leading 

asymptotics derived in (21 [20] is still valid. This has been proved in [13] by replacing the Witten- 
Laplacian by a supersymmetric Dirichlet operator and again using the Helffer-Sjostrand formula. 
In [31 [3] the authors were able to drop the assumption of translation invariance in the study of 
the correlation asymptotics. (At the same time they were able to relax the conditions on the 
on-site potentials of [21 [23 [21]-) They showed that, still for very small temperatures and for 
energy functionals with only one global minimum, the correlation asymptotics are determined by 
the Hessian of the energy functional at the global minimum. They derived this result for finite 
spin systems on graphs. Since all their results are again uniform in the number of spins, it is, 
however, easy to see that their methods work also in the infinite-dimensional case. 

In this article we choose to work in the thermodynamic limit and our goal is to analyse the 
formal Hessian of the energy functional evaluated at the global minimum more precisely. The 
latter is an example of a transition operator in So the main result of this article is actually 

to develop a method to calculate the asymptotics of the Green kernel of a certain class of operators 
in (^^). Instead of considering the large distance asymptotics we fix two points in R'^ and assume 
that they sit on lattices with vanishing lattice spacing, = {hx : x G Z'^}, h g (0, 1], d ^ 2. At 
the same time we assume that the spatial variation of the interaction potentials is fixed, too. In 
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this way we obtain a kind of continuum limit which permits to interprete the Hessian of our energy 
functional as a /i-pseudodifFerential operator on R''. We remark, however, that we introduce the 
continuum limit only for methodological reasons and that, in the limit, one would end up with 
a trivial theory. In order to analyse the Green kernel of the Hessian we therefore construct a 
parametrix for the /i-pseudodifferential operator. This is essentially done via a geometric optics 
-or WKB- construction for the solution of a corresponding heat equation. Here we are lead to work 
with complex- valued symbols by the introduction of exponential weights. Consequently, we shall 
deal with Fourier integral operators with complex- valued phase functions [111 [151 [TBI [19] . Since the 
results of [HI \W\ do not directly apply to our situation we show in detail how they can be suitably 
modified. In the end one verifies that the values of the Green kernel are, up to error terms of infinite 
order in h, given by the Fourier coefficients of the parametrix. Applying the method of stationary 
phase to these coefficients we find the leading asymptotics, as /i ^ 0, of the Green kernel. We 
obtain a formula which is similar to the expression known from the translation invariant case and 
which is determined by a certain Finsler structure naturally associated to the /i-pseudodifferential 
operator. If we consider the correlation, for two points in R'^ which are very close to each other, 
then we essentially observe again an Ornstein-Zernike behaviour. We remark that our results are 
closely related to those obtained in (HI [IE] where the authors develop a symbol calculus on rescaled 
lattices and discuss Finsler structures associated to pseudodifferential operators in order to study 
the tunnel effect for a general class of difference operators. 

The article is organized as follows. In Section [5] we state our hypotheses and our main re- 
sult precisely. In Section [3] we determine the /i-pseudodifferential operator mentioned above and 
conjugate it with exponential weights. In Section [3] we construct the appropriate exponential 
weight functions. In Section [S] we discuss the eikonal, or, time-dependent Hamilton- Jacobi equa- 
tion proceeding along the lines of [151 116] . The transport equations are considered in Section [6l 
In Section [7] we construct the inverse of the formal Hessian of our energy functional at the global 
minimum and, finally, in Section [5] we calculate the leading asymptotics. The text is followed by 
three appendices: In Appendix 1X1 we prove that the results of 3] are applicable in our situation. 
In Appendix [B] we present some elementary calculations which are used to reformulate the for- 
mula for the leading asymptotics in the translation invariant case. Some basic facts about almost 
analytic extensions are collected in Appendix ICl 



2 Hypotheses and main results 

For every h G (0, 1] and some d S N, d ^ 2, we consider a classical system of real-valued spins 
on the lattice := {hx : a; e Z'*}, which is endowed with the Euclidean distance. A spin 
configuration is given by an element of R^'' and denoted by cr = {'^x)xez'l- Furthermore, we write 
(jx = (o'x)xeT, for subsets T C Z^. Our lattice spin model is determined by the choice of a 
heuristic energy functional of the system. The latter describes the interaction of the spins and is 
(formally) given by 

0<\x-y\ih R 

Here, the coupling constant J > will eventually be assumed to be sufficiently small. R G [l,oo) 
is some arbitrarily large but finite interaction radius. We assume that the on-site potentials, D, 
and the pair interactions, W, fulfill the following 



Hypothesis. 2.1 (i) It holds D e C°°(IR'^ x R,R) and there exist C,C , x > and L > 2 such 
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that, for all {x,9) e R'^ x R, 

\d^D{x,9)\ ^ J/ = 0,1, 2, (1) 

dlDix,9) ^ -C, (2) 

\d^sDix,9) ~ d^gD{x,0)\ ^ C\dgD{x,9)\, (3) 

deD{x,9)sgn{9) ^ C m&x{\9\^+'' , \9\} . (4) 

For 9^0, we have d^D{-, 0) G Cj^iR'^) and 

2 inf d^D{x,0) > snp d'^D{x,0) . (5) 



(li) For allie , < |^| < R, we have W{-, (.,■■) G C°°{R'^ x R,R) and, for {x,9) eR'^xR, 

W{x,t,9) = W{x,-t,9) = ^/(a;,^,-^) , (6) 

W{x,l,9) ^ 0, Ty(a;,^,0) = 0, (7) 

|£| = 1 ^ a^W^(a;,£,0) > 1, (8) 

sup { |9^W^(x,£,6»)| + \dlW{x,l,9)\] < oo. (9) 
(2^,e)eiR''+i 

Kna%, dlW{-,l,Q) G C;;"(R''). 

Remark. 2.2 The estimates ([2p-f3i '^'"^ required in order to apply the results of T]\ \1S^. We 
assume that dgD{-,0) and dgW{-,£,0) are elements of C^iRf^), because this will allow us to 
interprete the formal Hessian of E at zero as a h-pseudodifjerential operator. We notice that 
Huvothesis \2.1\ implies that, for every x G Z^, the on-site potential D{x,-) has only one local 
minimum, namely a non- degenerate one at zero, and tends to infinity, as 9 ^ ±cxi, faster than 
any quadratic polynomial. We also point out that D{x, •) does not have to be convex. and ^ 
imply that d^W{x,£,0) ^ 0, for all x G R'' and £<eZ'',0 < \£\ < R. 



Example. 2.3 We let m G N, m ^ 2, and suppose that there are functions Cq, Ci, . . . , C2m € 
Cj^{R'^) such that inf C2m > 0, 2inf C2 > supc2, 

D{x, 9) = C2m{x) + ... + ci{x)9 + cq{x) , x G R'' , 61 G R , 

and deD{x,9) = ■i=^ 9 = 0, for x G R''. Furthermore, we suppose that W{-,£, ■■) = 0, for \£\ ^ 1, 
and that there is some w G C^{R'^), w ^ 1/2, such that 

W{x,£,9) = w{x)9^ , 

for \£\ = I, X G R"^, and 9 £ R. Then Huvothesis \2.1\ is fulfilled. If co,...,C2m and w are 

constant, we obtain examples of translation invariant models, which arise in lattice approximations 
of Euclidean quantum field theories of polynomial type and are known as "P(if)d-models" . 



Next, we recall the notion of a tempered Gibbs measure for the spin system under consideration. 
To begin with we introduce, for every finite subset T C Z^ , a (well-defined) local energy functional 
with boundary condition u G R^>^ , 



0<\j:-y\^h R 
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We further introduce a local specification, which is a family of stochastic kernels indexed by 
all finite subsets T C and depending on the parameters h £ (0, f ] and the so-called inverse 
temperature (3 € [l,oo). Given a tempered boundary condition, uj G UaeN^-Q' where — {a E 
R^i : ((x)" ax)x e ^li'^h)}^ their values are defined by 

f^l3.rji{A,uj) := I — ■ / 1a(o't|'^tO exp f- /3£^t(o't|'^tO') "^o-T , 



17'^ \ 



for A G S(IR^''), where Q5(IR^'') denotes the Borel-cr-algebra defined by the product topology on 
R^''. The partition function, w, ft.), is a normalization factor. For w G R^'' \ (UaeN^-a)' 

we set ^.p,^: ji{A,Lo) = 0. A probability measure, ^/3./t, on the Borel sets of R^'' is called a Gibbs 
measure (for the model determined by E) , iff it satisfies the Dobrushin-Lanford-Ruelle equilibrium 
equations, 

^lp^h,T{A,uJ) dfif3,hiuj) = Hf3,h{A) , T c finite, A G ^(R^') . 



A Gibbs measure fj.p,h is called tempered iff fJ-i3,h{f^'^a) — for some a > d/2. 

Under Hypothesis 12 . II the results of [T] are apphcable and ensure, for all /3 G [1, oo), J > 0, and 
h G (0, 1], the existence of a convex set, J, h), of tempered Gibbs measures. This convex set 
contains its extreme points, which are called pure, tempered Gibbs measures. We fix such a pure, 
tempered Gibbs measure, ^ip^h, in the following, let p denote the projection onto the constant 
functions in J^f^ :— L'^{^i3^fi), and set p-'- = 1 — p. We are interested in the asymptotic behaviour, 
as h 0, oi the truncated two-point correlation functions. 



C0Tf3.h{(Tx 
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for fixed x,y £ Z^, at very large inverse temperatures f3. In this situation we expect the spin con- 
figurations to be localized with high probability near the constant configuration {0)^^j^d, because 
all on-site potentials D(x, •) have a unique minimum at zero. In fact, it turns out that the asymp- 
totic behaviour of Corp h is determined by the formal Hessian of E at zero, E"{0) G ^(^^(ZJ^)), 
where E" {a) is given by 

{E"{a))^^ := d„^d„^E{y}{ay\a{yy) = {E"{a))^^, (10) 

for x,y € Zf and a G R^''. Here the second equality follows from which we also use to derive 
the following explicit formulas, 

{Uh{x), x = y, 

0<\x-y\^hR, (11) 
0, \x^y\>hR, 

for x,y £ J.'l, where 

Uh{x) D'Ux,0) + J W;>'g{x + ^,e,0), (12) 

o<\e\^R 

V{xJ) J<,(x,^,0), (13) 

for X G R'', < |£| ^ R. We state the precise relationship between £"'(0) and Corp,h in the 
following theorem which is essentially due to Bach and M0ller [4]. The latter authors actually 
consider finite lattice systems uniformly in the number of spins. Using the estimates obtained in 
[T5] it is, however, easy to see that their results extend to our situation. 

We remark that, due to (O and ©, we know that E"{0) is continuously invertible, if J > is 
sufficiently small. 
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Theorem. 2.4 Assume that Hypothesis \2.1\ is fulfilled and let ^ip^h G '^{P, J, h) he pure. Then 
there exist /3o ^ 1, Jq > 0, and C > 0, such that, for all (3 5^ Po, J € [0, Jo], and h G (0, 1], 

Proof: The claim follows from the results obtained in [5] together with the estimates derived 
in [IH] . We remark that all statements and estimates of these articles when applied to our situation 
hold uniformly in h since all constants in Hypothesis 12 . II are /i-independent. We present the details 
in Appendix [XI □ 

So, our aim will be to calculate the leading asymptotics, as /i — > 0, of the matrix elements 
(i?"(0)~^)_^^, for fixed x,y & Z^. The main point here is that we do not assume the spin system 
to be invariant with respect to translations on the lattice. In the translation invariant case, that 
is, if D and W do not depend on the first variable x, the asymptotic behaviour of the two-point 
correlation functions is well-known from Ornstein-Zernike theory; see, e.g. [T7]. We recall the 
precise result in Theorem 12.51 In the translation invariant case we can actually set h equal to one 
by a scaling transformation and consider the large distance asymptotics instead. The detailed form 
of the leading asymptotics (fTH|) given in Theorem 12 . 51 has been derived in [SlI^DlUlJ for translation 
invariant systems on finite discrete tori. As a starting point the authors used the Helffer-Sjostrand 
formula for the correlation which involves a certain Witten-Laplacian associated to the lattice 
model. The validity of (|16p for pure, translation invariant tempered Gibbs measures has been 
proven in [T^ again using a Helffer-Sjostrand formula. In the infinite-dimensional setting of the 
latter paper the Witten-Laplacian is replaced by a certain supersymmetric Dirichlet operator. 

Before we state the next theorems we have to introduce some notation. First, we define a 
Hamilton function, H e C°°(T*R'', R), by 

H{x,p) -.^ J2 V{x,e)e-<'^\P^ - Uix), (14) 

for a; e R'', p e Rd, where V is given by ((T^ and 

U{x) := D'^g{x,0) + J2 ^(^'^)- (15) 

For each x e R'\ H{x, •) is strictly convex and even and H{x,0) < 0, provided J > is small 
enough. It therefore makes sense to define the polar body 

K {peRd\H{x,p) ^0}, 

which is a strictly convex set that is symmetric about the origin and has a smooth boundary. 
Moreover, we let 

F{x,v) -.^ snp{{p\v)\peK} . ^'e^,R'* = R^ 
denote the support function of Ik*. For w ^ 0, we have the explicit formula 

F{x,v) = {p{x,v) \ v) , 

where p{x, w) e (Fj, := 91k* is the unique point, where the exterior normal field on (F^r points in the 
direction of v. In particular, F is smooth on the slit tangent bundle TR'^ = TR"^ \ (R*^ x {0}) and 
absolutely homogenous of degree one in v, i.e. F{x,Xv) = \X\F{x,v), for A e R, {x,v) e TR''. We 
also introduce the positive definite d x d-matrix G{x, v) with entries 

G(x, !))„■ := F{x, v) FX. {x, v) + F'^, {x, v) {x, v) , (x, v) e TR^ . 

Gij is homogenous of degree zero in u, i.e. G^j{x,\v) = G^j{x,v), for A e R, (a;, w) e TR"^. Finally, 
we let Hpp{x,p)^ denote the Hessian of H{x, •) at p ^ restricted to the orthogonal complement 
of VpJf(x,p) in Rd- In the translation invariant setting of the next theorem the functions H,F,p, 
and G are x- independent, whence we drop the reference to x in the notation of its statement. 
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Theorem. 2.5 Assume that Hvpothesis \2. 1\ is fulfilled, that J > is sufficiently small, and assume 
additionally that D{x,9) and W{x,l,6) are constant in x. Then, for sufficiently large (3, there 
is a unique tempered Gibbs measure determined by D and W. It is invariant with respect to 
translations on the lattice and its truncated two-point correlation functions fulfill 

^ f ^ 1 |V,gb(x-,))|"^ l + g(^)+0(r^/^) 



X exp ( - (1 + 
1 ^deiG{x-y) l + 0(^) + 0(r^/2) 



(3 {p{x~y)\VpH{p{x-y))) {2-kF{^))^ 
xexp(-(l + 0(/3-i/2))F(^)), (17) 

as (3 ^ oo and \x ~ y\/h oo, x,y ^ if^. 

Proof: Using (jl6p . which is derived in the above mentioned papers, and some elementary 
hnear algebra we deduce Formula (fTTl) in Appendix [B] □ 

In the following theorem we state the main result of this article. To formulate it, we introduce 
the Finslcr distance associated to F, 



dFix,y) ^ M J F{q,q) , x,y&R'^, 



where the infimum is taken over all piecewise smooth curves q : [0, t] R'* such that r > 0, q{Q) — 
y, q{T) = X. It turns out that the infimum is always attained by at least one (smooth) geodesic from 
y to X. Since F is absolutely homogenous of degree one in v any reparamctrization of a minimizing 
geodesic yields again a minimizing geodesic and dp is symmetric. We can always reparametrize 
any minimizing geodesic in such a way that we obtain the projection of a Hamiltonian trajectory 
defined by H which runs in the level surface (F := -ff~^({0}). If there is, up to reparamctrization, a 
unique minimizing geodesic from y to a; we denote its initial (resp. end) velocity and momentum 
in the just mentioned special parametrization by Vy and Py (resp. and p^). We remark that 
we have {px \ vx) > and Vx = V pH{x,px) and analogous statements for y. As in Riemannian 
geometry one may define the notion of a conjugate point for x and y with respect to F; see 
Section [3] below. 

Finally, we set /i„ := 1/2", n e N. If TV e N and x e Zf^^, then of course x G Zf^, for aU 
integers n N. 



Theorem. 2.6 Assume that Hypothesis is fulfilled and that J > is sufficiently small. Let 
N E x,y G J.'l , and assume that, up to reparamctrization, there exists a unique geodesic 
minimizing Finslerian arc length from y to x and that x and y are not conjugate to each other. 
(This assumption is always fulfilled, for fixed y, provided x is sufficiently close to y.) Then 



^det {G{x,Vx)G{y,Vy)) exp {- dF{x,y)/h„) 
A{x,y) y/{px\vx){py\vy) {2t, dF{x,y)/hn)^ 

d+l 

+ 0{hn^ ) exp{ - dF{x,y)/hn) , 



as n oo, where 

A{x,y) = 1 + Oidpix^yy/^), as x ^ y . 
Proof: The theorem is proved by combining Propositions 14. 5r ii). 17. 11 [Ol fc 18.21 and Equations 

(iiniii&dnsi). □ 
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Remark. 2.7 The function A{x,y) can be expressed explicitely in terms of transversal Jacobi 
fields along the unit speed geodesic, q, from y to x; see Proposition 1 8. 2\ and the remarks preceeding 
it. If the geodesies emanating from y are dispersing along q, then A(x, y) increases when dpix, y) 
gets large. If they are bunching together, /S.{x^y) decreases and vanishes, if x is conjugate to y. If 
the flag curvature \SI happens to be zero along q, then A{x, y) is equal to one. 

3 A related /i-pseudodifFerential operator 

The basic idea underlying our analysis is the fact that, for every / e £'^{Zf^), the expression (recall 
(HBl-dni)) 

{E"{0)f){x) ^Ur,{x)f{x) - (18) 
can be interpreted as the image of any / e C°°(R'^) with f\-id= f under a /i-pseudodifferential 

h 

operator evaluated at the point x G Z^. Before we explain this correspondence in detail, we fix 
some notation and recall some general facts. 

For our purposes it is sufficient to consider only bounded symbols: Let f2 be an open subset 
of R'^, R'^ X Rd, or [0,oo) x R'' x etc., where T'' := {R/2tiTY is the d-dimensional torus. Then 
we write / G Sb{n), for a function / : x (0, 1] ^ C, iff /(• ; h) e C°°{n), for aU h € (0, 1], and 
there is some Hq G (0, 1] such that, for all multi-indices a, 

sup \d^f{uj;h)\ < oo. 

(w,/i)Gax(0,fto] 

We recall that the Weyl quantization, Op)^(a), of any symbol a G «S't,(IR'* x Rd) is determined by 
the oscillatory integrals 

Opr(«)/»= /e^<«l^-^>/''«(^,e)/(j;) feyiR'), xGR^ (19) 

where S^{R'^) is the space of Schwartz test functions on R'^. The possibility to regard (|18p as the 
action of a pseudodifferential operator is a consequence of the following: If a G >S'&(IR'' x Rd) is a 
symbol which is (27rZ)'^-periodic in ^, then its Weyl quantization has the distribution kernel 

Kaix,y) = (^-ia)(^,x-2/) = /i'^ ^ a(^, z) (5,(y - x) , (20) 

where u){S,) — J e^*^^ ' ^^Z'* u{x) dx is the semiclassical Fourier transform and 

aix,z) = ^^e-<^l«)/V^,0 zeZi, (21) 

are the Fourier coefficients of a{x, ■). So if / G ■y{R'^) we get, 

OpJ!'{a)f{x) ^ J2 a{^,y-x)f{y), x(.R'. (22) 

yGs+Z^ 

To apply these remarks to E"{Q) we consider the Fourier series corresponding to the V {■,€), 

V{x,0-^ E e*^^"^na:,^) = V{x,l)co^{{(,\l)), xGR^eeCrf. (23) 

0<|f|^-R 0<|«|!S-R 

(Here the second expression follows from the first since V is even with respect to t.) As a trigono- 
metric polynomial V{x, ■) is of course entire and (27rZ)''-periodic with respect to real ^. By 
Hypothesis 12. 11 [//,,, V{-,1) G ^^(IR'^), and we see that the symbol Uh — V is contained in the class 
Sb{R'' X Rd). 
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Proposition. 3.1 For all f e ^^(ZfJ and f e C°°(R'*) with /tzd^ /, 

{E"{0) /) {x) = Op^iU, - V) fix) , X e Z^ (24) 

Proof: Since {V)''{x,z) = h-'^V{x,z/h) = h-'^V{x,-z/h), < |z| < hR, and = 
otherwise, p4p is a special case of (\22^ with x € Z^, provided / S oS^(IR''). We know, however, 
that Op^ {Uh — V) maps C°°{R'^) into itself, because only finitely many Fourier coefficients of V 
are non- vanishing. Consequently, ([22|) is available, for every / G C°°(IR'^), in this case. □ 
Proposition 13.11 suggests to study the inverse of E" (0) by means of the semiclassical microlocal 
analysis of the operator Op^ {Ut — V). To work out the exponential decay of {E"{0)~^)xy, we 
shall, however, first conjugate the operators in fM]) with suitable exponential weights and then 
construct a parametrix for the conjugated operators. Therefore, we assume that we are given some 
weight function (p G 5';,(IR'^) in the following and use the symbol e'^^'^ also to denote the diagonal 
multiplication operator on i?^(Z^), whose action on the canonical orthonormal basis vectors is 
given by e'^^'^Cx ■= e'^'^'^^^'^ex, x € Zf. We remark that, eventually, it will be necessary to consider 
also conjugations with exponential weight functions that depend additionally on a time and a 
momentum parameter. 

In the next lemma we derive a classical asymptoic expansion of the Weyl symbol corresponding 
to the conjugated operators. As we shall see in the proof this is actually possible using only 
completely elementary arguments since we are dealing with trigonometric polynomials. 

Lemma. 3.2 There is a classical Weyl symbol, a^ E >5'b(IR'* x IR^j), such that 

e"^/"- Opr([/,. - V) e'^l^ fix) = I Opr(at) f{^) , ^ e , (25) 
and, in particular, 

e^/" E"{0) e-*/'^ fix) = t Opr(a^) f{x) , xelf,, (26) 

for all f G C°°(IR'') and f G ^^(Z^) such that /tz^= /• '^tv ^ trigonometric polynomial in ^ 
and admits a classical asymptotic expansion in 5*6 (R'* x IR^), 

oc 

aU^S ^Y.h'^'aU^^O, (27) 

i/=0 

where a^^ai^, . . . are trigonometric polynomials in too. Its principal symbol is given by (recall 

m) 

a'I'ix, e) := 4ix, = i Vi^, ^ + i(t)'ix)) - i U{x) , {x, ^ e R'^ x . (28) 

If (j> = (j>{t, X, rf) G C°°([0, oo) X R'' X Rj) is such that (j) + i{x\ri) defines an element of Sb{[0, oo) x 
ipd jd-j^ ^^g^ symbols a^, 00,03,... can be viewed as elements of Sb{[0, oo)t x R!^ x IR^^j x T^) 
and the expansion (27\ ) is valid in Sb{[0,oo) x R'' x x T''). 

Proof: First, we consider the multiplication operator e"^ ^ Op^ (Uh) e~'^/'^ = Uh- Using ©, 
which implies V{x,£) = V{x, —£), for all x G R'^ and < |£| ^ R, we write 

Unix) = DUx,0) + I J2 {Vix + f,i) + Vix~fj)). 

0<\l\^R 

By Taylor's formula, we find, for every e N, some '^^{■^i^h) G Sh{R'^) such that 

V{x + ^,i) + V{x-^J) = ^-_(-) \i\^)^-V{x,i) + h^''+^-^N{x4-h). 
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By (I1I1)-(I131) we further have, for all x G R"^, 
e"^/^ Opf (y)e~'^/'' f{x) 



Here the integrals under the double sum are non-zero only \i y ~ x — hi. If the latter condition 
is fulfilled we have, however, again by Taylor's formula, 

\ {m - m - I (<^(^ - ¥) - ^(^ + ¥)) (30) 

= V0(^) ) - — ^ (I)'' ( ^ I ^ r+'cj^C-^) + h^^^^ $^(£±H, ^ ; /,) , 



^ (2J/+1)! V2 

where $Ar(-,£; K) G Sb{R'^), for every G N. Plugging (|5D|) into writing out the exponential 
series 

N ,2 

( - E (^) < ^ I ^ )^''+ V(^) + <i>A.(^ , ^ ; /.)) (31) 

and ordering the terms with respect to powers of h we thus get, for every € [Kl, 



e^/''Opr(f/,.-y)e-^/V>) (32) 



where Oq , , ■ ■ ■ and r2, r4, . . . are Weyl symbols, which are trigonometric polynomials in ^. For 
the principal symbol we find 

atix,0^^ E V{xJ)e'^'\^^-^'\'^'^-^'> - zU{x), 



which is (UHl)- Applying ^ and ((211) to ^ we obtain (f^ -d^f l) . 

If now (j) = X, T]) is such that (/) + i{x\r]) is 27r-periodic in 7y, we observe first that dtcj), (j)'^, 
and d^^ ^ ^^(f), for |a| ^ 2, are again 27r-periodic functions. For the first derivative with respect to 

X we get (j>'^{t,x,r] + e) = 4>'^{t,x,r]) — ie, for e S (27rZ)''. This shows that all derivatives of 
appearing in (PTjl are 27r- periodic in as well as e-<^ I for £ G Z'*, < |^| ^ i?. Therefore, 

all symbols Oq , 03 , • ■ • and r2, r4, . . . can be viewed as elements of 5'f,([0, 00) x R'* x IRj; x T''), if 
this is possible for (f> + i{x\ri), and it is clear that the last statement holds true. □ 

In order to construct a parametrix for Op)f (a^) we will solve a corresponding heat equation 
by means of a WKB construction. For this purpose we state the following lemma, which is a 
simple consequence of the method of stationary phase. We suppose again that (f) depends on two 
additional parameters such that (j) + i{x\ri) G S'6([0, 00) x R'' x IR^) is 27r-periodic in 77. We further 
introduce the vector field 

d 

Y{t,x,ri) :=^%a*(*-'')(x,C)|^^oa,., 

and set af = af = ■ ■ ■ = 0. 

Lemma. 3.3 Assume that b^, G S'f,([0,oo) x R'' x T"*), for v G Nq; and suppose that there is some 
compact subset, JC d R"*, such that supp {bu{t, •, 77)) d JC, for G Nq and {t, rj) G [0, 00) x T'^. Let 
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b(t, x,rj]h) he a Borel resummation of the formal series X^^^o h'^b^(t, x, rf). Then we find symbols 
Fat e 5'6([0,oo) X R'^ X T''), e N, with 

supp {rN{t,-,ri)) <Z K, + {\x\ R) (t, 77) e [0, oo) x , 

such that, for all {t,x,ri) e [0,oo) x R'* x T'' and sufficiently small h> 0, 

(zOpr(at)6)(i,a;,r;) 
= m'^(*'-''')(a;,0)6(t,a;,r;;/i) + h {(Y + ^divY)bo){t,x,r]) 

N-l v-2 

+ $]/i"(((i^+^divr)fe,_i)(t,a;,77) + 6^)(t,a;,77)) 
+ h'^ riq{t, x,r];h), 
where, for f & C°°(R^), 

4 Construction of a weight function 

In order to study the exponential decay of the Green kernel of E"{0) we distinguish two points 
X* and y* and design a suitable weight function, if, which models the exponential decay of 
(E"{0)~^) ^ ^. Here we can choose any pair of points x*,y* satisfying a certain condition in- 
troduced below. An appropriate weight function will essentially be given as a solution of the 
Hamilton- Jacobi equation 

T{x,Lp'{x)) - U{x) = , (34) 

where T{x,p) := V{x,ip). For our purposes it is actually sufficient to solve (|34[) only in a small 
neighbourhood of a certain Finsler geodesic from y* to a;*. The solution ip will then be extended 
in such a way that, for all x outside that neighbourhood, it holds T{x, f'{x)) — U{x) < 0. The 
construction of ip is presented in Subsection 14.11 below. In Subsection 14.21 we collect various 
properties of the principal symbol a"^ defined by with ip = ip, which play an important role 
in the sequel. 

4.1 The construction of if 

There is no reason to restrict ourselves to the special Hamilton function T — U appearing in l|34p 
for the purpose of this subsection. In Subsection 14.21 we will observe that T — t/ is an element of 
the set of all Hamilton functions, H, that satisfy 

Hypothesis. 4.1 It holds H £ C°°{T*R'^, R). For all x E R'^, the function H{x, ■ ) : T^R'^ R is 
strictly convex, even, and H{x,0) < 0. 

We thus seek for solutions of the Hamilton-Jacobi equation H{x,ip'{x)) — 0, where H satisfies 
Hypothesis 14.11 The appropriate solution will essentially be given by the Finsler distance function 
to some prescribed point on the prolongation of the geodesic from y* to x* . Before we turn to the 
construction of ip we recall some required notions and facts from Finsler geometry. In Section [5] 
we shall also make use of them. All the standard results collected below can be found, e.g., in 

[siaiii]. 

We have already introduced the polar bodies 

B<; {peRd = T;R'^ : H{p, x) } , x e R'^ . 
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For every x £ R'^, Ik* is strictly convex and symmetric about the origin. Its boundary, 

{peT^R" : Hip,x) - 0} , 

which is cahed the figuratrix at x, is a smooth submanifold. We further set f :— H^^{{0}). We 
also recall that the support function, F{x, ■) : R'^ — T^R"^ — > [0,oo), of Ik* is given by 

F{x,v) -.^ sup{{p\v) : pekl} , v e R"^ . (35) 

The set 

K := {veR'^^T^R'' : F{x,v) < l} 
is again strictly convex and symmetric about the origin, and its smooth boundary, 

:= {f eT^r' : F{x,v) = 1} , 

is called the indicatrix at a; £ R'^. For ^ 0, the value of the support function is given by 
F(x,v) = {p{x,v) I u), where p{x,v) is the unique point on at which the exterior normal field 
on tx, i.e. the normalized gradient of H with respect to p, equals v := v/\v\. The function 
F is continuous on TIR'* and smooth on TR'^. (Here and in the following we write TR'^ for the 
tangent bundle with the zero section deleted and use a similar notation for the cotangent bundle.) 
Moreover, F{x,-) is absolutely homogenous of degree one, i.e. F{x,0v) — \9\F{x,v), for all 
{x,v) E TR'^ and 9 E R. Euler's theorem thus implies 

Fl{x,v) = p{x,v) and F^'^{x,v)v , {x,v)EfR'^. (36) 

The restriction of F!^^{x, v) to the orthogonal complement of v is strictly positive. The matrix G 
with entries 

is therefore positive definite at every point of TIR''. Altogether F enjoys all properties required 
for a Finsler structure on TR'^; see, e.g., [5]. In the calculus of variations F is called an elliptic 
parametric Lagrangean and H is an associated Hamiltonian in the sense of Caratheodory, as 
H'p^Oon t*R'^ and H{x, F^{x, v)) = 0, for (a;, v) E TR'*; see, e.g., [B]. We also introduce 

Hypothesis. 4.2 

inf { F{x, v)\xER'^ , vE S"*"^ } > . 



Since F is absolutely homogenous of degree one in v, the map dp '■ R"^ ^ R'^ ^ [0,cxd), given by 

dF{x,y) := M^iq), J^{q) := J F{q,q) , (37) 

defines a metric on R'^. Here the infimum is taken over all piecewise smooth curves in R"^ from y 
to X. (One could equally well take the minimum over Lipschitz continuous curves defined on the 
unit interval, which would give the same results.) In [111 §5.3] it is shown that 

lp{x) :^ dFix,y), x E R"^ , (38) 

is the (unique) viscosity solution of 

H{x, ip'{x)) = , for X e R'^ \ {y} , ^{y) = . 

In particular, this means that tp is continuous, differentiable on R** \ ({y} U E), where S is some 
closed set of finite {d — l)-dimensional Hausdorff measure [TU], and satisfies H{x,ip' {x)) = 0, for 
X ER'^\ ({y} U S). Moreover, it is known that, for any function u : R'^ ^ R with u{y) = 0, which 
is locally Lipschitz continuous and satisfies H(x,u'{x)) ^ almost everywhere, it holds u ^ip. In 
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Section E] it will become clear that the latter property shows that ^ is an optimal weight fmiction 
for our purposes. Since it is not smooth we have, however, to modify it suitably. To this end we 
recall some further notions from Finsler geometry in the following. 

Extremals of the functional are called geodesies. It is well-known that under Hvpothesis l4.2l 
the infimum in ((57)) is actually attained and that each geodesic is smooth. The homogenity of F 
implies that si is invariant under reparametrizations. In particular, every reparametrization of 
a geodesic is again a geodesic. Moreover, we have the following well-known lemma relating the 
geodesies to trajectories of the Hamiltonian vector field 

Its proof can be found, e.g., in pp. 197]. 

Lemma. 4.3 Any geodesic defined by means of F can be reparametrized in such a way, that we 
obtain the projection onto R'' of a Hamiltonian trajectory defined by means of H and running in 
(F C r*IR'^. Conversely, any projection of a Hamiltonian trajectory running in ^ is a geodesic. 

Next, we recall that under Hvpothesis l4.2l the exponential map determined by the Finsler structure 
F is defined on TR'' by exp^(u) := q^A'^), for [x, v) e TIR'', u 7^ 0, and exp^(O) := x, for x g R''. 
Here qx,v is the unique geodesic which passes through x at time i = with velocity v. (It holds 
Qx.viOt) = Qx.evit), 9 > 0.) If w G Oa;, then 1 1-^ qx,v{t) = exp^(i w) is the unit speed geodesic passing 
through a; at t = in the direction of v, that is, F(^qx,v{t),qx.v{t)) = 1, for all t £ R. We remark 
that the projections of the Hamiltonian trajectories onto R'^ might have different velocities. 

We assume that x,y,v G R'', and x = expy{v) in the following. If the derivative of exp^ : 
TyR"^ R"^ is singular at v, then x is called a conjugate point for y. If x is the first conjugate 
point for y along the geodesic t i— s- q(t) := expy{tv), then we know that, for all r e (0, 1), (ztfcr] 
is, again up to reparametrization, the unique minimizing geodesic among all geodesies from y to 
X that run in some small neighbourhood of q{[0,T]). It does not have, however, to be globally 
minimizing. For r > 1, ^[[cr] will definitely lose its unique minimizing property even among nearby 
geodesies. The cut point of y along q is by definition q{r), where r is the supremum of all those 
T > such that q\[o^T] is globally minimizing. The cut point of y along q always appears before 
or at the first conjugate point. Notice that, if x is simultanously the cut and first conjugate point 
for y along q, then (ztici] might still be globally minimizing. The unique minimizing property of q 
will, however, definitely fail at the cut point, if the latter occurs strictly before the first conjugate 
point. 

We can always find some open set, W C TyR"^, which is star-shaped with respect to zero, such 
that expj^f^-e C'^{W)nC°°{W\{0}) is bijective from W onto := expy{W), and it is well-known 
that 

ip{x) = dF{x,y) = F{y,{expy\y^)-'^{x)) , x e ; 
see, e.g., Chapter 8, §3.3]. Moreover, we have 

^'{x) = F^(y,7(r)) = ^(r) , x&'^\{y}, (39) 

where t ^—f {l{t),S,{t)) S (F is the Hamiltonian trajectory corresponding to the minimizing geodesic 
from y to X such that 7(0) = y and 7(t) = x, for some r > 0. 

In the following we fix two distinguished points x*,y* E R'', and assume 

Hypothesis. 4.4 x* ^ y* . Up to reparametrization, there is a unique minimizing geodesic from 
y* to X* . X* and y* are not conjugate to each other. 

We recall that Hypothesis 14.41 is always fulfilled provided x* and y* are sufficiently close to each 
other. Since F is absolutely homogenous there exists of course, up to reparametrization, again 
only one minimizing geodesic from x* to y* if HvDOthesis l4.4l is satisfied. Moreover, we can prolong 
the geodesic from y* to x* or from x* to y* a little bit such that it is still minimizing. 
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Proposition. 4.5 Assume that H e C°°(T*R'*, R) fulfills HypotheseslJl\ &\T^ and that x* , y* G 
R'* fulfill Hypothesis \4-4\ Then there exist a point, yg, on the prolongation of the geodesic from 
X* to y* , a compact neighbourhood, Kq, of the geodesic segment from y* to x* , some open set, 
W C TygR'', which is star-shaped with respect to zero, and some function, ip G C^(R'',R), such 
that the following holds true: 

(i) For all x £ R"^, 

H{x,p'{x)) «C 0, 

H{x,ip'{x)) = xeKo. 

(a) (p{x) — (p{y*) — dp{x,y*), for all x on the geodesic segment from y* to x* . 
(ill) exp^J^re C^{W) n C°°{W \ {0}) is mjectwe on W and 

C exp^„(#')\{yo}. 

(iv) For every x G Kq, there is a unique pair (r, po) G (0, oo) x such that the projection of 
[0, r] 3 t i—> exp{t XH){yQ,Po) onto R'' is a minimizing geodesic from yo to x. We have 

exp{TXH){yo,Pa) = {x,ip'{x)). 

Proof: Let / be an open interval and g : / ^ R"* be a geodesic with q{ti) — y* and (7(^2) = x* , 
for some ti,t2 G /, <i < <2- We can find some to G /, to < ^i, such that q\[to,t2] s^^^^, up to 
reparametrizations, the unique geodesic segment from yo '■= q{to) to x* . We can moreover find 
some open set W C Ty^R'^, star-shaped with respect to zero, such that the restriction of exp^^^ to 
W is bijective from W onto := expy^(W), where g([io,i2]) C . We have exp^^ G C^^) H 
C°°{W \ {0}). Defining ^ hy §EI with y = yo we know that ip{x) = F{yo, (exp^^^ 
X e-W, solves H{x, ip'ix)) = on \ {yo}- 

Next, we pick some open set, and two compact sets, Ki, and K, with 

q{[to,h]) ^ C <£ K <Z K ^ C <£ . 

We suppose that the boundaries ^K^, and dK are smooth submanifolds. Furthermore, we can 
assume that exp~^^(i4r^), cxpy^{K), and expy^^(^j,) are star-shaped with respect to zero. We pick 

some H G (R'', R) with ^ = on K^,, and < H{x) -H{x, 0), for a; G R'' \ if*, and set 

:= H + e H , e G (0, 1] . (40) 

For every e G (0, 1], we let (p^ denote the function constructed exactly in the same way as ip using 
the Hamiltonian Hf, instead. We notice that = ^ on K^,, since the Finsler structures and, 
hence, the geodesies emanating from yo corresponding to and H are identical on Kt,. 

In the following we show that, for all sufficiently small e G (0, 1], ip^ is, apart from yo, smooth at 
every point in a neighbourhood of K. We know that all geodesies of the Hamiltonian H contained 
in and emanating from yo form a central field about yo. No conjugate or cut points of yo are 
present in . First we claim 

Lemma. 4.6 There is some £0 € (0, 1] such that, for all e G (0,So], every geodesic of the Hamil- 
tonian emanating from yo leaves Ki, and touches a point of R"* \ '^i, before it possibly enters 
Ki, again. 

Proof: Let t; G 8^^,, and let r > be the first exit time of Ki, for the geodesic 1 1— > expy^(t u) 
defined by H. Then the geodesies t ^ ex\)y^{tv) and t 1— exp^^(iw) are identical on [0, r], for 
£ G (0,1]. Here we designate the exponential map defined by by a superscript e. We also 
know that t ^ exp (tv) touches a point in R'* \ say, at time t > 0, before it possibly enters 
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Ki, again. There is an open neighbourhood, J^^ C ly^^ of v in the indicatrix at such that 
e-yi^y^itw) e ^'^\'^ , for aU w e Since all geodesies depend smoothly on the parameter e and 
on ui e J^i,, we find some e„ € (0, 1] such that eyi])y^{Jw) € R"* \ for all w G J^i, and e G (0, £i,]. 
The claim now follows from the compactness of ly^. For the open cover contains a finite 

subcover IJ^^j^ J'v, and we can set £o min{e.„j, . . . , e^^}. □ 

Lemma. 4.7 There is some ei > such that, for all e G (0, ei], there is no conjugate or cut point 
for yo contained in K , if the geodesies are defined using the Hamiltonian H^. 

Proof: The assertion about the conjugate points is clear since they are characterized by 
singularities of the exponential map, which depends smoothly on e. Concerning the cut points 
we prove a slightly stronger statement, namely that there are ei G (0, 1] and c > such that, 
for all e G (0,£i], and every pair of unit speed geodesies, qf and g|, qf ^ g|, defined by and 
emanating from yo, we have 

(3t,i>0 : ql{[0,t]) C X A ql{t) ^ ql{i)^ =^ t > i + c. 

Note that, for unit speed geodesies, the time parameter equals the Finsler arc length. We argue 
by contradiction and suppose that there are sequences, (£„), (c„), (i„), (t„), (w„), («„), such that 
En \ 0, c„ \ 0, as n — > oo, and > 0, Vn,Vn € '^x* , Vn 7^ Vn, for n G N, which satisfy 

exp^;; ( t„ «„) = exp^;; ( in v„) , and exp^;; ( s G X , s G [0, i„] , 

for n G N. By definition of K^, we know that the unit speed geodesic with initial velocity Vn 
does not stay inside /-iT^, for all t G [0,t„]. By Lemma [4.61 we thus find r„ G (0,t„) such that 
exp^jj(T„ Vn) G for ri G N. By compactness and by a choice of suitable subsequences we can 
assume that the sequences (in), (wn), (vn), and (t„), have a limit, as n tends to infinity, 

which we denote by t,t,v,v, and r, respectively. Then it holds t ^t, eKpy^{Tv) G '^Z, 

exPj^o(*^) = *=^Ps/o ( ^~'^) ^ ^^^^ expj^^(sw) G , s G [0, i] . 

It also follows that v v and we get a contradiction, because [0,t] 3 s ^—^ expy^(sv) is the unique 
minimizing geodesic for H from yo to exp^^^ ^) • 

Returning to the proof of Proposition HTSl we notice that Lemma H77I implies that, for all e G (0, £1], 
is smooth in a (in general £-dependcnt) neighbourhood of _ftr \ {yo}j since the set where (pe is 
not smooth is closed. 

Next, we pick some compact set, K', satisfying Ki, K' C K' d K, and some cut-off function 
X G C°°(IR'', [0, 1]) such that x = 1 on /C' and x = on K". We note that there is some C > 
such that — H ^ C e on supp (x')- We further pick a family, {j,5}5>o, of Friedrichs moUifiers 
such that supp (^,5) (e {|a;| ^ (5} and observe that 

{xVe + {l-x)is*'Pe)' = XVe + {'^ - X) 3S * "f'e + ("Pe - js * ^e) x' (41) 

holds on R'^ \ {yo}- Here we note that the expression [js * '■p'i;){x) makes sense, also for x G K'^, 
since tp^, is differentiable almost everywhere. Since the function H appearing in the definition (|40p 
of iJg depends only on x, the set 

K {pe Rd : H,{x,p) 0} 

is still convex, and by the properties of H it holds lk| C Ika;, for all x G R''. This inclusion is 
proper if a; G K'^. Moreover, we have H^{x,ip'^{x)) = 0, that is, (p^(a;) G lk|, for almost every 
X G R''. Since js is a probability density and k% is convex, it holds [js * ^'e){x) G lk|, for all x G R'*. 
Consequently, for every x G R''\{yo} the convex combination x{x) v'six) + {'^~x{x)) {js*(p'g){x) is 
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contained in IkJ, too. On the support of x' the set lk| and the indicatrix have a strictly positive 
distance which is uniformly bounded from below by some contant C{e) > 0. We also know that 
the term ((^^ — js * (fe) converges to zero uniformly on the compact set supp(x'), as 5 tends to 
zero. Choosing S > small enough, we can therefore ensure that the right side of ([1T|) is contained 
in Ik:^, for all x g IR'*\ {yo}- 

The function x fe + ~ x) js{e) * Ve is smooth everywhere except at uq. To obtain a smooth 
function we pick some Q G C°°(IR, [0, oo)) such that ^ O' ^ 1, 8 = const > in a neighbourhood 
of zero, Q{t) — t, for all dp{yo,y*)/2 < t < 2dF{yo,x*), and 9 = const, on [3dp{yo,x*),oo). 
Finally, we choose some e G (0, £i], and some sufficiently small 6{e) > 0, and set 

e o ( X + (1 - X) js{e) ■ 

We observe that we still have (p'{x) G Ik^;, for x G R**, since ^ 0' ^ 1. To complete the proof we 
note that, by the minimizing property of the geodesic from yo to x*, 

ip{x) - (p{y*) = dF{x,yo) - dpiy* ,yo) = dF{x,y*) , 

for all X on the unique geodesic from y* to x*. Finally, (iii) is fulfilled by construction and (iv) 
follows from □ 



4.2 The principal symbol a'^ 

We turn to the discussion of the special Hamilton function H := T — U , where 

T{x,p) = V{x,ip) = ^ cosh((£|p)) , xeR'^,peRd, (42) 

and U is given by (jl5p . First, we observe as in [20l §4] that, for every a; G R'', the function 
T{x, ■) : Rd (0,oo) is strictly convex and even. Moreover, < T{x,0) ^ U{x) — g, for some x- 
independent g > 0, provided J > is sufficiently small. In other words, H satisfies Hypothesis STTJ 
Using the uniform bounds on U{x) and V{£,x) of Hvpothesis 12.11 it is not difficult to see that 
Hypothesis 14.21 is fulfilled by H, too. So, Proposition 14.51 is applicable, for each pair of points 
x*,y* satisfying Hvpothesis 14.41 

Since we want to use the results of Section [3] with (j) = ip, we collect some important properties 
of the principal symbol 

a'P{x,C) ^ iV{x,i + iip'{x)) - iU{x), (x, C) G R'' x R^ , (43) 

in the next lemma. In its statement and henceforth we use the short-hand notation (a'^)^'^. '■— 
dx^ ^a^ etc. The Hamiltonian vector field of Sfta'^ is denoted by Xs^q.^ := (.^^j^J^)- 

Lemma. 4.8 Fix two distinguished points x*,y* G R** fulfilling Hypothesis \4-4\ o,i^d let and Kq 
be as in Proposition \4.5\ Then the symbol a"^ G S'fc(R'' x R^;) has the following properties: 



(i) For {x,£,) G R'' x R^, 

^a'^ix,^) 0, (44) 

5a^(a;,^) = ^ (x,^) e Kq x {2TrZf . (45) 

(li) For (a;,0 e Kq x (27rZ)^ 

{a'^Kx,^ = 0, (46) 

{a%{x,0 = (5fta^)^(x,e) - if;(x,^'(a;)) ^ 0, (47) 

[a^'A^^O = 0, (48) 

{an'li{x,0 = -iH';A^,p'{x)), (49) 

{a'^rU^^O = H';Ax,^'ix)) + H';A^,p'{x))if"ix), (50) 

= 0. (51) 
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(Hi) p : [0,r] — > Kq x (27rZ)'^ wii/i t > is an integral curve of Xsrq*., i/ and only if, for 
{x{t),£,{t)) := p{t), the curve [0, r] 3 t i-^ (^x{t), tp' {x{t))^ is a piece of an integral curve of Xh 
emanating from some point on {yo} x (Fj,o C r*IR''. 

(iv) Let 7 : [0,t] Kq denote the geodesic from y* ~ 7(0) to x* ~ "/(t) parametrized such that it 
is the projection onto R'^ of a trajectory of Xh ■ Then, for every constant vector e £ (27rZ)'^, (7, e) 
is an integral curve of Xs^av ■ 

Proof: (i): Similarly to analogous considerations in we deduce (|33]) and (|45p from the 
representation 

V{x,C + ip'{x)) = ^ y(rE,£){cosh((£|^'(.T))) cos((£|0) 

0<|^|^i? 

-zsinh((^|^'(x))) sin ((^ 10)}, (52) 

the positivity of V, and the inequality H{x,tp'{x)) — V(x,i(p' {x)) — U{x) ^ 0, a; G R'^. In the 
derivation of ([^5)1 we make use of the assumption that V{x,£) ^ 1, if \£\ — 1, and the trivial fact 
that, for all ^ & Rd \ i'^nZ)'^, we find some k e Z'^, = 1, such that cos((fc | C)) < 1- 

(ii) follows from straightforward computations. In particular, we obtain, for all x G Kq and 
e G (2^Z)^ 

{na^y^{x,0 - E V{x,fYismh{{i\^'{x))) = T;{x,p'{x)) + 0, (53) 

0<|f|<H 

To see that the terms in ((55)) are in fact non-zero, for (x,^ G ifo x (27rZ)'', we can refer to the 
strict convexity of T(x, •), or multiply ([55)1 with (/^'(a;) and use again that F(x, ^) > 0, for \i\ — 1. 

(in) Let p — (x,^) : [0,t] Kq x {2'kTY be a trajectory of X^af- Then (|T7)) implies 
i;(t) = V5a'^(p(t)) = Vpi?(x(i),(/?'(a:(t))), for aU t G [0,r]. Moreover, differentiating the identity 
H{x,ip'{x)) — 0, X e Kq, we get ip"{x)VpH{x,ip'{x)) — —V xH{x,Lp' [x)), x G Kq. From this 
we infer ■^if'{x{t)) — ~VxH{x{t),ip'{x{t))), and we see that [0,t] 3 t t-^ {x{t), tp' {x{t))) is an 
integral curve of Xh- By Proposition 14.51 we know that, for every x G Kq, {x,ip'{x)) lies on an 
integral curve of Xh emanating from some point on {yo} x ^y^. By the unique solvability of 
Hamilton's equations we conclude that [0, t] 3 t 1-^ ^'i^it})) is a piece of an integral curve 

of Xh emanating from {yo} x [Fj,^. 

Conversely, if [0, t] 3 1 1-^ (2^(^)7 is a- piece of an integral curve of Xh emanating from 

{yo} X fyo! then we use ()47p and (|15|) to check that [0,r] 3 t 1-^ (^x{t),e) is an integral curve of 
XsRa, for every e G {2nZ)'^. 

(iv) follows from (iii) and Proposition 14. 5r iv). □ 

5 The Hamilton-Jacobi equation 

In this section we take a first step in the construction of a parametrix for Op)j'^(a^), where if is 
the weight function constructed in Proposition and is given by Lemma 13.21 Since we shall 
obtain the parametrix by integrating a parametrix for the corresponding heat operator we seek 
for approximate solutions of the time-dependent Hamilton-Jacobi equation 

dtij{t,x,7^) + a'^{x,i>',{t,x,f]}) ^ 0, i;{0,x,7j)^{x\T]), (54) 

where t G [0, 00), x G R'^, and rj G Rrf. Here and a'^ is given by ([28| with (j) = ip. Since the symbol 
is complex-valued we have to work with almost analytic extensions and we can only hope to 
solve (IMl) up to error terms of order ON{\'^'ip{t, x, ?7)|^), for every G N. Approximate solutions 
of have been constructed in 15J for a class of symbols which arc homogenous of degree one 
and in |16j for purely imaginary symbols satisfying some further assumptions. Since neither of 
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these works applies directly to our situation we will explain in detail how the proofs of 15j can 
be adapted to fit with our hypotheses. We remark that the same class of problems is adressed in 
[3] where approximate solutions to the heat equation are obtained in the framework of Maslov's 
canonical operator theory. Approximate solutions of complex eikonal and transport equations are 
also derived in [^ . 

This section is split in two parts: In Subsection 15.11 we extend the crucial estimates from jl5j 
on the real Hamiltonian flow associated with to our situation. The essential properties of the 
symbol required here are that the imaginary part of a"^ has a fixed sign and that the derivative 
of ^a"^ is non-zero if Q'a'^ vanishes. Both of them are ensured by our special choice of f as we 
have already observed in Lemma l4.8l In Subsection 15 . 21 we construct the approximate solution ij). 
Some basic facts about almost analytic extensions are collected in Appendix O 

5.1 Estimates for the Hamiltonian and contact flows 

We will obtain the approximate solution, ^/;, of (I54p essentially as the generating function of the 
flow of a real Hamiltonian vector field associated with . In order to control the error terms 
and to assure that the imaginary part of if) has the right sign it is necessary to derive suitable 
estimates on the flow of this vector field. This is carried out in the present subsection by adapting 
the arguments of [E], where the symbol is assumed to be homogenous of degree one. The same 
problem has also been considered in [9]. We remark that the reasoning of [15] and, hence, of this 
subsection works globally, too. 

In the following we denote the complex coordinates in T*£,'^ = x Cd again by (a;,^) = p, 
so that d^i = ^{dsn^i — id^^i), = ^(SsRa;; + idc^^i), etc. We pick almost analytic extensions of 
(/?, [/, and V < |£| < i?, and denote them again by the same symbols. Plugging them into 
P5)) we obtain an extension of a"^ to C'' x Cd which is almost analytic with respect to the first and 
analytic with respect to the second variable. More precisely, a''^ is a trigonometric polynomial in 
^ with almost analytic coefficients. We denote this extension of a'^ simply by a in the rest 
of Section [5l 

Since a as well as any partial derivative of a satisfies the Cauchy-Riemann differential equations 
on the real domain we know that a'p[p) = a'yf^p{p), a'pp{p) = '^ytpUp(p)^ ^'^^ every real p, and, hence, 
(|T7ll -(fFT |l still hold for {x, S,) £ KqX (27rZ)'^, if the derivatives are interpreted as complex derivatives. 
We emphasize that (|44| and ((45|) are available only for real (a;,'^)- 

We are going to study the fiow of the real field [1^] 

ATa := Xa + Xa , 

where 

d 

1=1 

We verify that 

d 

Xa = E{5R(4J%.. + - - SKO^^eJ, (55) 

1=1 

so that the definition of Xa amounts to considering Xa as a vector field on R^''. (If we took 
another almost analytic extension of a then the components of the corresponding real fields would 
be equivalent in the sense explained in Appendix [C]) We denote the flow of Xa by 

Ktip) {Qit,p),E{t,p)) exp(tXa)(p) , peT*C''. 

On account of ([55)1 we see that the Hamiltonian equations for Q and S read 

Q^V^a{Q,E), S = -V,a(g,S). (56) 
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Due to the periodicity of a in 3?^ we have, for {t, y, 77) e IR x T*C'^ and e G {2TiTf, 

Q{t,y,V + e) ^ Q{t,y,v) , '^{t,y,V + e) ^ E{t,y,'q) + e . (57) 

It will be convenient to extend T*C'^ by an extra variable, s, which parametrizes the action. We 
may view the resulting space C x T*C'* as a contact manifold equipped with the contact form 
ds + X^iLi dx"^- Then we may call the vector field 

Ka Ka + Ka, 

where 

Ka := -Ads + Xa, A -.^ {^\a'^) - a, 

the real contact field associated with a. Here and in the following ( • | • • ) denotes the bilinear 
extension to of the Euclidean scalar product on R''. The function A is called the elementary 
action. To study the flow of Ka we introduce the function 6 : C x T*C'^ — > C, where 

6{s,x,0 := -3s - (5x|3fiO , {s,x,0 e C x T*C'^ . 

The relationship between these objects and the ones considered in [TS], where the symbol is 
assumed to be homogenous of degree one, can be easily seen by means of the the standard reduction 
from non-homogenous to homogenous symbols: Denoting the momentum variable conjugate to 
the extra variable s by a, we set a'^°™(s, x,a,£,) := a a{x, ^ C)) fo^' (^7 ^7 '^^ ^ T*^^^^ with a ^ 0. 
The results of [15] apply directly to the flow of the real field X^hom , which is defined analogously to 
Xa- Since a is obviously constant along the flow lines of X^hom it is easily checked that X^jhom t{(T=i} 
can be identified with the real contact field Ka- Moreover, the function — ( 3(s, x) \ 3?((j, ^) ) used 
in [15j reduces to S, for cr = 1. (A direct application of [15j would, however, give suboptimal 
versions of the estimates (|5^ -(fM |) below.) 

Lemma. 5.1 For all compact subsets f2 d T*^*, there is some Cq G (0, 00), such that, for all 

{s,p) eCxn, 

Ka{&){s,p) ^ -l^a{dip) - Cnl^pf. (58) 

Proof: Taylor expanding at dip in the third and fourth step, we get 

Kai^s + (3a;|3?C)) 

= 3(a(p) - (^|4(p)) + (e-*3^l4(p))) - (5a;|3?<(p)) 
= 3a(p) - (5Ra;,(p)|3p) 

= 3a(5Rp) + i (3a;:,(5Rp)5p| 3p) + 0(13^^ 
= ^5a(3fip) + i (3a(3fip + 2 3p) + Qa(5Rp-2 5p)) + 0(13^1^) 
<; ^ cja(5Rp) + 0{\^p\^) . (59) 

□ 

Remark. 5.2 Exactly as in '15', §1], we may check that 6 is invariantly defined modulo 0(\'^{x, ^)|^) 
on compact subsets with respect to coordinate changes x '—^ x. (Of course, the covariant variables 
^ have to be transformed accordingly.) Then the arguments of U5\ Remark 1.8] directly imply that 
Estimate {5<$l) does not depend on the particular coordinate system. 

In the following we fix some po = (a^OiCo) S T*IR'^. We assume that there exists some complex 
neighbourhood, C C^, of xq and an almost analytic function, -00 G C°°{y) such that ^ipo ^ 
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on y n R'^, SV'oCa^o) = 0, and iP'q{xo) = S,q. By a calculation similar to (l5^ wc then see that, for 
every real, compact subset d n R'^, wc find some compact y 'f" and some C £ (0, oo), 
such that y'nR'' = and 

3Vo(a;) - (9a;|3fiV'o(a;)) ^ -C'|3x|^ x € . (60) 

We introduce the set 

£o {(-V'o(a;),a;,V^^,(a;)) : a: G r } , 

which we may call an almost analytic Legendre submanifold. (In our applications below we will 
only consider the case where ^'0(2^) ■— "0(01 2;, 77) := {x\r]), for some real rj e Rd-) Furthermore, 
we assume that there is some r ^ such that nt{po) = exp(tXa)(po) is real, for all t £ [0, r]. 
If r > 0, we additionally assume that Sa(po) = 0. In this case it holds (3a)'(Kt(yOo)) = 0, and, 
hence, ^a{Kt{po)) = 0, for t £ [0, r]. By we also know that (5Ra)'(Kt(po)) 7^ 0, for t £ [0, r]. 
Denoting the flow of Ka by 

$4 exp{tKa) 

we have that 

*t(s,p) = (Q(s,p),Kt(p)) , q(s,p) := s - / A(Kr(p))dr. (61) 

Jo 

In particular, — ■*/'o(a;o), Po) is also real, for t £ [0, r]. 
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Lemma. 5.3 In the situation described above there exist e > 0, C £ (0,oo), and some neighbour- 
hood ff <Z £, y. T*C'^ of {—tpo{xo) , Po) such that the following inequalities hold on £0 H ^, for all 

!^r s^t + e, 

1 rt 

(62) 

(63) 
(64) 

In particular, if (s,p) g £9 H ^ and Kt{p) and q(s,p) are both real, for some t £ [0,t + e], then 
Kr{p) is real for all r £ [0,t]. 

Proof: Taylor expanding the right side of ^ Qnt = QXai^t) at SKkj and using Duhamel's 
formula we obtain, exactly as in [T^ pp. 351], the estimate 

ISk^I < 0(l)(|9Kt| + f \^a'{nKr)\dr'] . (65) 

It holds for all p in a compact, complex neighbourhood of po and u,t £ [0, t + Ei], for some ei > 0. 
Since, for p in a compact set, the curves [0,t + Ei] 3 t ^Kt{p) stay in a compact set, we may 
apply the standard estimate for positive functions to (|65p . which together with Holder's inequality 
gives 

I'^Krl 0{l)(^\'^Kt\ + (y -'^a{^Ky)dvy^^Y (66) 

for ^ M < r ^ t ^ T + £1. Next, we integrate the estimate ([SH) for JjS(<i>f) = Ka{&){^t) from 
u to t, use (j66p to bound r £ [u,t], and arrive at 

6($t) > 6i^u) - \ f '^a{Kr)dr - 0{1) {t - u) {'^ml^ - ( /" -"^ai^Kr) drY^"^ . (67) 
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For u = 0, (pU]) implies 

S($o(-V'o(x),a:>^(x))) ^ 

for x in a sufficiently small neighbourhood of xq . Here the last term can again be estimated by 
means of ([TO]) . We further notice that we can make the integral appearing in ((S7| arbitrarily small 
by assuming that p is contained in a sufficiently small neighbourhood of po, and that £i > is 
sufficiently small. For, ii t G [0,t], the integral vanishes at po. Setting u = we thus get 
The proofs of ([63]) and ([64]) are identical to those in [151 PP- 355], except for Equation (3.30) in 
[15], where — 9q has to be added on the right side. □ 

5.2 Approximate solution of the Hamilton-Jacobi equation 

The aim of this subsection is to construct an approximate solution, ijj, of ([5^ . We will define 
■0 by a formula well-known from classical mechanics and proceed along traditional arguments to 
check that ■0 has the required properties. In doing so we use Lemma 15.31 to controU the error 
terms and to ensure that 3'0 is positive. Since Lemma 15.31 is applicable near real integral curves 
of Xa, we start by considering the latter more closely. All the time we keep on using the notation 
of Proposition [i31 with H = T-U. 

Lemma. 5.4 Let r > and p : [0,t] T*R'^ be a real integral curve of Xa with p(Q) E {3a = 0}. 
Then p{[0,t]) C Kq x (27rZ)'*. If we write p{t) = {x{t) , S^(t)) , then [0, r] 9 (x{t) , ip' [x(t))) is a 
piece of an integral curve of Xh emanating from some point on {j/o} x fyo T'*IR'^. Conversely, 
every point of Kq lies on the projection of an integral curve of Xh emanating from some point on 
{2/0} ^ fyo '^'^'^ piece of that projection inside Kq is the projection of a real integral curve of 
Xa with momentum e E {2ttJ.)'^. 

Proof: p being a real integral curve of Xa means (pit)) = ^ ^pit) = 0. Since a fuffiUs 

the Cauchy-Riemann differential equations on the real domain, it follows that (5a)!ffp(p(i)) = 
'^{a'p){p{t)) — 0. Hence, the derivative of (9a) |"t*ir'' vanishes along p and 3a(p(0)) = thus 
implies p([0,r]) C {3a = 0}. Using (liS)) we conclude p([0, r]) C Kq x {2'kTY. The remaining 
assertions follows from Proposition 14.51 and Lemma [THT iii'l. □ 

Of course, there might be real integral curves of Xa defined on a non-trivial interval running in a 
region where 3a < 0. These do, however, not play any important role since the imaginary part of 
the action is strictly decreasing along them. 

The approximate solution ^ will turn out to be "almost" a generating function of the canonical 
relations 

Ct = {{x,^,y,i^)ET*€^xT*£^ : (x, = «t(2/, ^) } , t^Q. 

We recall our notation = exp(iXa). For every t ^ 0, the relation Ct is 27r-periodic in the sense 
that (x, ^, y, 77) E Ct implies (a;, ^ -f e, y, 77 + e) e Cf , for every e E (27rZ)'*. To study Ct we denote 
the bilinear extension of the canonical symplectic form on T*IR'' to T*C'^ by ct, i.e. 

ct(z9, = *^ J 0^)^'' (^TT*C'^ . 

Furthermore, we write 

CTeCT(z9,^') ;?((79,,^^), «,z9^)) - ct((79„z?,), (^;,<,)), 
where § = {§^,i)^,^y,i)rj) e T{T*C'^) x T(T*C'*) and similarly for i?'. 

Lemma. 5.5 Let t > and p : [0,r] — > A'o x (27rZ)'' be a real integral curve of Xa- Then Ct is 
strictly positive along p, that is, for all t E [0,t], 

l5eCT(7?,^) ^ 0, 7?eT(p(,),p(o))Ct, (68) 
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with equality if and only if -d is contained in the complexification of the tangent space at {p(t),p{0)) 
of Ct n R''. The latter condition holds, if and only if = 0, dy, 0). 

Proof: Along the real curve p{t) =: {x{t),e) G Ko x (27rZ)'' the fundamental matrix of a is 
given by 

where 

A{t) H';j,{x{t),^'{x{t))) , Bit) := H';,{x{t), ^' {x{t))) + A{t) ^" {x{t)) , (70) 

for t e [o,t]. It holds ^ = ie{t),eio)) = {e^{t),e^it),e,io),e^io)), where e{t) = {e^{t),e^{t)) - 

k;(p(0))6'(0), t e [0,r], for some ^(O) e T^(o)r*C'', and we know that 0{t) satisfies f^9{t) = 
fa{p{t)) 0{t), for t e [0, r]. We thus obtain 

I 1 5(0(0,0(1)) = (3?0,.(t)|A(t)5R0at)) + (3%(i)|A(t)30aO), 

for i e [0,t]. Integrating the previous identity from to t and using the fact that A(i) is positive 
definite, for all t € [0,t], we see that assertion of the lemma holds true. In fact, the structure of 
fa{x{t), e) shows that we have equality in (|68p if and only if %(0) = 0. □ 

In the following we set 

^ { (t,y,e) e [0,c5o) X ifo X (27rZ)'' : k*. (y, e) G i^o x {e} , t' G [0, t] } , 

& := {{t,x,e) e [0,oo) X Ka X {2TrZf : K^t'{x,e) e Kq x {e} , t' e [0,t]} , 
S ({0} X R'' X Rd) U ^ , := ({0} x R'^ x R^) U ^ . 

By Proposition 14.51 and Lemma [5.41 it is clear that there is some t^ e (0,oo) such that QJ^'S Q. 
[0,to] X X (27rZ)'' and, moreover, 

V (t, y, e) e ^ : Q(t, y, e) = ^pH{Q{t, y, e), (^'(g(i, y, e))) . 

Lemma. 5.6 There is some neighbourhood, 'i§ , of S in [0,oo) x x with the following prop- 
erties: There are two smooth functions, k and g, defined on 'S , such that, for every {t, x, £,,y,ri) g 
[0,oo) X T*C'' X T*^^, 

{x,^,y,v)eCt 44> y ^ k{t,x,ri) A ^^g{t,x,r]). 

is l-K -periodic in dirj and, for all {t, x,rj) E '^^ and e £ {2ttJ.)'^, 

g{t, a:, 77 + e) = g{t, x, 77) + e , k{t, x,ri + e) = k{t, x, 77) . (71) 

Moreover, 

V(i,a;, e)G^: g(t,x,e) = e, k{t,x,e) = tt e^p {~ t Xh){x, (p' (x)) . (72) 

Proof: It suffices to show that the projection Ct 3 {x,£,,y,ri) 1-^ {x,rj) is bijective in a 
neighbourhood of each point on Ct with (t, x, rj) G For t = 0, this is trivial by definition of Cq. 
So, it remains to treat the case {t, x,ri) € t > 0. 

Let t > and let p : [0,t] Kq x {2ttZ)'^ be a real integral curve of Xa- Let = 
{'dx,'&i,'dy,'d7j) € T(^p{t),p{o))Ct be such that dx = and = 0. It clearly holds a Qa{d,d) = 0. 
Using the same notation as in the proof of Lemma [5751 we thus have that %(r) = 0, for all r G [0, t]. 
In particular, i)^ ~ 0. Then it is also clear that = 0, because 6x{r) solves the initial value prob- 
lem -§fOx{r) = IB(r)6'^(r), Oxit) = 0. Consequently, the projection Ct 9 {x,£,,y,ri) 1-^ (^,77) is 
bijective in a neighbourhood of {p{t),p{0)) G Ct- 
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Finally, ([7T|) follows from the periodicity of Ct and ([7^ from Lemma [53] and the definition of 



In the remaining part of this section we show that an approximate solution of ()54p is given by 

^{t,x,T]) := {k{t,x,T])\7j) + f A{Kr{Ht,x,rj),Tj))dr, {t,x,Tj)e'^. (73) 

Jo 

We observe that ^(;{0,x,r]) = {x\i]) and a direct computation using ([SS]) . ([ST]) . and ([7T|) shows, 
for {t,x,'n) e ^ and e e (27rZ)'^, 

■0(i,a;,77 + e) = V(i,a;,r7) + (?7|e) , (74) 

so that the difference -0 — (a; | r/) is (27rZ)'^-periodic in 5Rry. By definition of and ((45|l and (|47p 
we also know that 

V (i, x, ry) e : !3V(t, a;, ??) = . (75) 
To study tjj further we introduce, for {t,y,r]) G IR x r*C'^, 

z{t,y,r]) := -Qi-{y\v),y,v) ^ {y\v) + / A{i^riy,v))dr , 

Ja 

f{t,y,r^) := |3(g, S)(t, y, ry)]' + 6(-Z, Q, S)(t, y, r?) . 

Lemma. 5.7 There is a small neighbourhood, jV , of S in [0, cx)) x x IR^, which is (27rZ)''- 
periodic in the last variable, such that, for all {t,y,Tf) g and r G [0, i], 

( 3Q(t, y, 77) 1 3?S(i, y,rj)) - ^ 3a(5RK„(y, 77)) (76) 

-0(l)|3(Q,S)(t,y,r;)|', 
r(r,y,77) < 0(l)f(t,y,77), (77) 

where the O-symbols are uniform on compact subsets of .jV . 

Proof: For each fixed 77 G IR^, we apply Lemma [5.31 with 7/'o(j/) := (yl??): y & C"^, and 
po e X {v}- If Po e i^o X (27rZ)'*, we choose r := max{t > : Kt(po) G i^o x (27rZ)''}. If 
Po ^ -K^o X (27rZ)'', we set t — 0. So, if t > 0, then Lemma [5.41 ensures that Kt{po) is real, 
for t G [0,r], and by ((45|) we know that 3a(po) = 0. Then ([76|) follows from ([62]). The first 
inequality in ([77)1 follows from ([^ . ([75]) . and the fact that we can assume |3((3, S)(r, y, 77) | to be 

arbitrarily small by restricting our attention to some sufficiently small neighbourhood of where 
S)(r, y, 7/) I vanishes. The second inequality in ([77]) follows from ([U]) and ([^5]) . □ 

In the next lemma we compare the differential of Z with the pull-back under of the Cartan 
form, 

uj := ^dx — a{x,^)dt, 

considered as a form on IR x T*C'^. Following a standard proof known from classical mechanics 
(see, e.g., pp. 479]) we only have to keep track of the error terms. 

Lemma. 5.8 (i) On every compact subset ^ IR x T*C'^ such that \t\ ^ to on f2, 

dZ = -e*Lu + ydr] + Oa, ( max f (i, y, 7?)^) . (78) 

\iKto 



^Z{t,y,'q) ^ 
S(Q,S)(r,y,77)|' < 
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(ii) Let jV he as in Lemma \5.7\ Then 

dZ = -Q*uj + ydrj + iC'Ar(r^) 

on ,y/ , where the Oj^ -symbols are uniform on compact subsets of -VV . 
Proof: We set \:= dZ + 'd*uj and use 

z = (s|4(g,s)) - a(g,s) = (s|q) - a(g,s) (79) 

to obtain 

A = (Z;-(S|Q;))dy + (4-(S|Q^))#+ (Z;-(S|Q;))d,7 + {ZL^-{^\Q'^))d^. 

Now, let a denote any of the variables y'- ,rii,r]^, 1 ^ i ^ d, and set := {Z'^ — 

Using successively ([7^ . the Hamiltonian equations fS^ and the almost analyticity of a we find 

A„ = d^{{E\Q) - aiQ,E)) - {E\Q'J - {E\Q'J 
= (S'J4(Q,S)) -d4a{Q,0) + {aUQ,^)\Q'a) 

= 0^(|5(Q,S)r). 
Due to the initial conditions {Z, Q, E)\t=o — (( 2/ 1 ?7 ), ?/, '7) we have 

Z'y\t=o — V 1 Z'rj\t=o — y , Z^lt^Q — Z^lt^o — , Q'y\t=o = ^, Q',j\t=o = Q^\t=o — Q^\t=a = . 

We conclude ^ 

X{t,y,r,) ^ ydr, + f S)(r, 77)!^) dr , 

Jo 

which implies (i). (ii) follows from Lemma [5.71 (Note that (|77p is available only for real 77.) □ 
In what follows we set 

K{t,x,'q) := {t,k{t,x,i]),r]) , 

for all {t, X, rf) G so that 

K*Z = ZoK = ij, K*Q = QoK = X, K*E = EoK = g. (80) 
Furthermore, we write := K~'^{,A') and 

r := K*f = ToK = \^{x,g)\'^ - (9a;|3?g) + V- 
The set c/K C is a neighbourhood of in [0, 00) x x IR^, which is 27r-periodic in 77. 

Proposition. 5.9 Let ip be given by i7S\ ). Then 

dit,..^,^r,){dt^ + <^ix,^',)) - 0^,„(r^), (81) 

dit,.,^.,^r,M - 9) = OnA^^), (82) 

dit,.,^,n^M - k) = Ojv,a(r^), (83) 

= Ojv,a(r^), (84) 

= Ojv,a(r^), (85) 



on for € N oTid every multi-index a G Nq''^"'^. On ,a-symhols are uniform on compact 
subsets of -jV . 
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Proof: On account of Lemma [5751 11) we obtain 

= d{K*Z) = K*{dZ) 
= K*{q*uj + ydr^ + tON{y^)'^ 
= gdx - a{x,g)dt + kdr^ + tON{T^) 
on .yV . Using Lemma [CTD we Infer that ([5T|l - ([55|l hold true. □ 

The next corollary summarizes the properties of ip on the real domain, where the weight F can 
actually be replaced by 9-0: 

Corollary. 5.10 (i) There is some neighbourhood, of S" in [0, oo) x IR"* x IR^;, which is 2tt- 

periodic in the last variable, such that, for all {t,x,ri) £ 

^Ht.^.v) ^ ^\^9{t.x,ri)W (86) 
'^'>l){t,x,ri) = <^ {t,x,ri)eS'. (87) 
(Here the O-symhol is again uniform on compact subsets of ^n.) Consequently, 

so that \81\) - W^) hold true on with the right sides replaced by OnH^iP)^)- In particular, 

d(t,.,w(^*'J' + = ON.Ai^^f) on . 

(a) For all (t, X, e) G 

dt4'{t,x,e) = 0, '^ri'>P{t,x,e) = n exp — t X h) {x , (p' {x)) , tp'^[t,x,e) = e. (90) 

Proof: (1) We observe that (gl]), ([7S]), and ^ Imply 

(3a;|5R.g) ^ -0{\)\'^{x,g)\^ on ^. 

Let (to, 2;o, 770) G S . From (j72p we know that vanishes on S and we can thus find some compact, 
real neighbourhood, K' , of (io, a;o, ryo) in .J>^ f] ([0, 00) x R'' x R^), such that x + eQg{t, x, rj) S =yf , 
for all {t, X, rf) G K' and every e g [0, 1]. Consequently, there exist constants C, C S (0, 00) such 
that, for all {t,x,rj) e K' , 

'i'ip{t, X — e^g{t, x,ri),r]) ^ —C\'ig(t,x — £^g(t,x,r]),r])\'^ ^ — C' |3g(t, x, 77)]'^ . 
Taylor expanding the left hand side with respect to x and using ([5^ and ([51|) we obtain 

for some A^o G N, A^o ^ 2, and C" ,Cmo & (0,oo). Now, we choose e S (0,2^) such that 
£ Catq I^VI^"^^ < 1/2 on iiT'. By possibly restricting our attention to an even smaller real compact 
neighbourhood, K" ^ K' , of (tg, 2:0, 770), we may further ensure that C"|9g| < e/4 on K". Then 
we readily obtain ([86|) on if". 

Concerning ([ST)) . It only remains to prove the Implication "=>" because of ([75)) . So, suppose 
that (i,a;,r7) € .^r and "^tpit, x,f]) = 0. Then ([511) Implies that Qg{t,x,r]) = and d??]) and PJ]) 
show that (y,??) := {k(t,x,r]),ri) and {x,g{t,x,ri)) are connected by a purely real Integral curve 
of Xq. Hence, If we had '^a{y,rf) < 0, then we also had Qtp{t,x,r]) > by ((ii)) and l(7^. So, we 
must have 'ia{y,r]) = and Lemma 15.41 Implies that i = or {t,x,ri) G 

(11) follows from (gB]), ([721), and (l81|). □ 

To solve the transport equations In the next section we need another consequence of Proposi- 
tion [ 
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Corollary. 5.11 Locally on JV , 

k\^^^ = OA,(r^) , g\^^^ = OA.(r^) , (91) 

and locally on some small neighbourhood, J*/' C , o] S in [0,oo) x x 

Q\w) = ^^(f"^) ' = ^^A^(r^) , (92) 

Proof: follows immcdiatly from ([5^ - ((551) . The first equality in ([^ is easily derived 

from ()9ip by differentiating the identity k{t,Q(t,y,r]),r]) = y with respect to y and {jj,rj). To 
obtain the second we differentiate g{t,Q(t,y,ri),r]) — E{t,y,r]) with respect to (y, fj) and use the 
first and again (PT|) . □ 



6 The transport equations 

In this section we carry out the next step in the construction of a parametrix for the heat operator 
hdt + Op^ {a^,). We seek for a parametrix, Q, of the form 

Qf{t,x) = J^^^'^^^^y'^-^<y^^^/H{t,x,7j;h)f{y)^^, (93) 

where tp is given by (|73p . For the amplitude h we use the ansatz 

oo 

b{t, X, 77 ; /i) X ^ h^b^it, x, r,) . (94) 

To find the sequence of transport equations that determines b(),bi, . . . we ignore for the moment 
that ip is defined only in a small neighbourhood of and argue formally. (We shall perform 
everything properly in Proposition 16.31 where the results of this and the previous section are put 
together.) Observing that 

e-^^/hOp^{a'^)e'^/'' = 6"^'^/''+^/'' Op)f (i F - i C/,,) e^^/''+'"^/'' = OpfAaV^) 
and using Lemma 13.31 with (j) — Lp — we see that the equation 

(/^a, + Opr(a^))e'^^/n = 
leads to the following sequence of transport equations: 

(To) (at + r + idivy)6o = Oiv(r'^), iVeN, 

(Ti) (5t+r+idivr)6i = P2bo + On{T''), N€U, 

v-l 

(T,) (5t + y + idivy)6, = + OA,(r^), TVeN, 



Here we have 



(95) 
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and the operators P2 , ^3, • ■ • are given by with (p ^ ip~ ii}^. By ([7i|) ip — iij} — i{x\ri) depends 
(27rZ)''-periodicaUy on 77 and, hence, the same holds for Y and so that 77 can be viewed as an 
element of T'^ in the following. As initial conditions we choose 

60(0, x, 77) := x{x), (96) 
6^(0, x, 77) 0, i/eN, (97) 

where rj €1'^ and x S C^(IR'', [0, 1]) is supported in some small neighbourhood of Kq and fulfills 
X = 1 in another small neighbourhood of Kq. In the following we also assume that x is extended 
almost analytically such that the extension, is compactly supported in some tiny complex 
neighbourhood of supp (x) ■ 

Lemma. 6.1 The transport equations (Tq), (Ti), . . . with initial conditions WOjl and |ff7| ) are solv- 
able in a small neighbourhood, JY' C jV , of S" in [0, 00) x x Rj_, which is 211 -periodic in 77. 
More generally, 

/or every € N and any multiindex a G Nq'^"'""'^, where the On .a- symbols are uniform on compact 
subsets of ,yV' . The solutions b^, v G Nq, are 2TT-periodic in 77. We have supp (5j/(i, •, 7;)) C 
Q(t, supp (x), 7;) , for all {t, 77) anrf G Nq- 

Proof: The proof is essentially standard; see, e.g., [551 PP- 578]. We have, however, to keep 
track of the error terms. We argue by induction and consider the equation 

(dt + Z)b ^ fb + w, b\t=o = c 

where 

d 

Z{t,x) := ^d^^a'^{x,g{t,x,r]))d^i . 
1=1 

Moreover, / := — ^ divK and w is either 0, in which case we set c = x, or w = X]^=o P^+i->t 
where we suppose that &01 • ■ • , ^i^-i have already been constructed. In the latter case c — 0. We 
know that 

/4 = O^(r^) on ^ 
and as an induction hypothesis we assume that 

(&o)^, . . . , [K-i)'^ = Oiv(r^) on ^' R-^Jy"') , 

which implies 

w'-^ = C'Ar(r^) on JY' . 

Here .yK' is the neighbourhood of S appearing in CoroUarv lS.llI We observe that Q is the flow of 
Z := Z + Z. We first solve the following initial value problem for an ordinary differential equation, 

dtb{t,y,ri) = f{t,Q{t,y,Tj),r])b{t,y,T]) + u;(i, Q(t, y, 7/), 77) , 6(0,?/, 7;) = c(y) , (99) 

and observe that 

b{t,x,'q) := b{t,k{t,x,i]),r]) 

fulfills 

(dtb + Zb){t,Q{t,y,ri),ri) = ^b{t,Q{t,y,7j),T]) = f{t,Q{t,y,r]),r])b{t,y,r]) + w{t,Q{t,y,r]),r]) 

(100) 
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on jY' because of /c(t, Q(t^ y, rj), rj) — y. Next, we consider the error term 

b'-{t,x,v) = by{t,k{t,x,ri),T])l4:{t,x,ri) + by{t, k{t, x,j]),t]) k'^{t, x,!]) . (101) 

Here the error can be controlled using (PT|) . To investigate by we differentiate (IM)) with respect 
to y and obtain 

dtbL(t,y,7^) - f{t,Q{t,y,ri),r,)bLit,y,r,) = Ojv(f (i, y, 77)^) on , (102) 

if we take into account. It also holds b^{0,y,r]) — c^{y) — OjvdSyl^), so we infer from ([77)1 
and l|102p together with Gronwall's lemma that 

bL{t,k{t,x,7^),Tj) = ON{T{t,x,rif) . 

In view of we obtain 

6^(t, X, 77) = Ojv(r(i, X, 77)^) , (t, X, 77) G ^' . 

Together with ()100|) this shows that 

9t6 + Z5 = /6 + u; + C'Ar(r^) on (103) 

Finally, we observe that the components of Y and Z satisfy — Zi = C'jv(r^), because of (15^ . 
The vector field Z may hence be replaced by Y in p03p . 

Finally, it is clear that b is 27r-periodic in 77. In fact, the maps /, w, Q, and k used to 
construct b only depend on the equivalence class of 77 in T**. Moreover, it is easy to see that the 
inclusion supp {w(t, •, 7;)) C Q{t, supp (x); ??) implies that the support of b{t, •, 7;) is again contained 
in Q{t, supp (x), 7?)- So, we also obtain the last statement of the lemma. □ 

Next, we consider bo{T, x, 0) more closely, for (r, x, 0) G By definition of ^ we know that there 
is some y £ Kq such that Q(t, y, 0) = x and Q{t, y, 0) G Kq, for all t G [0, r]. By (071), jSS]), plj) . 
and (|95p we further know that, for t G [0, r], 

y(t, Q{t, y, 0), 0) - VpH{Q{t, y, 0),ip'{Q{t, y, 0))) = Q(t, y, 0) . (104) 

In particular, Q{t, y, 0) and Y{t, Q{t, y, 0), 0) are real. 

Henceforth we denote the determinant of a real square matrix, M, by \M\. 

Lemma. 6.2 Assume that (r, x) G ^. Then the solution of (To) with initial condition ^96]) 
satisfies 

bo{T,x,0) ^ \Qy{T,y,0)\ x{y)i where y = tt^ exp(-r X//)(x, (^'(x)) . 

Proof: The proof is well-known. We recall it for the sake of completeness. From the proof of 
Lemma HI] and (fT(H|) we infer that, for (t, 7/,0) G ^, 

^feo(i,Q(t,tJ,0),0) = -i (6odivy)(t,Q(t,y,0),0). 
The Liouvillc formula shows that 

^ I Q'y{t, y, 0) \-"^ = -i I Q'y{t, y, 0) f'^' (div F) (t, Q(i, y, 0), 0) , 
and it follows that 

6o(t,Q(i,y,0),0) = |Q;(i,y,0)P'/'x(2/)- 

□ 
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In the next proposition we draw the main conclusion of the results attained in Sections [5] and [H] 
so far. 

Up to now the restrictions of the coefficients b^{t, x, rj), v d No, to the real domain are defined 
in some 27r-periodic neighbourhood, fl of S" which is contained in the domain of i/j- In 
what follows we consider them as functions on [0, oo) x R'' x T'' by simply setting 6^ equal to zero 
outside -yV' n We pick some smooth cut-off function, g G C^([0,oo) x R'' x T'^), such that 
p = 1 in a neighbourhood of S' and supp (g) C n .-#ir, and let b denote a Borel resummation 
in 5'b([0, oo) X R'' X T'*) of the formal series on the right side of 

oo 

b(t,x,r];h) x h" x, rj) g{t, x, rj) . 

By Proposition 16. II we then know that there is some > and some compact set Jiif <e such 
that supp (6) C [0, to] X ^ X T''. 

We may assume that we find another compactly supported cut-off function, g E C^([0, oo) x 
ipd ^ jd^ jg^ j^j-j^ such that ^ = 1 on supp (6) and supp (g) C Extending ip again by zero outside 
we define 

ilj{t, X, 7]) := g{t, X, ry) ip{t, x, J]) + (1 - g{t, a;, 77)) (( x 1 77 ) + i) , 

for t ^ 0, a; e R'^, and 77 £ Rj.. tp is a smooth function such that Qip > on ([0, 00) x R'^ x R^) \ (o 
and such that ip{t, a;, 77) — ( x 1 77 ) can be viewed as an element of S'6([0,oo) x R'' x T"*). 

Proposition. 6.3 For every TV g N, there is some compactly supported G <S'6([0, 00) x R"* x T"*) 
with supp {rN{t, •, rj ; /i)) C J(f + {l^;] ^ R} such that, for sufficiently small h > 0, 

{hdt + Opr(a^)) (e^^'/H) = h'^e'^/^fM. 

Proof: By definition of af^ with 4> = (p and (p = ip — itp, respectively, 

^~^^/h Opf (a^) e'^/'^ b = e-*^/'' Op)f(a^) e^^/'' b = Op^^ia"^'"^) b . (105) 

By our choice of ip the phase tp — iip + i{x\r]) defines an element of 5b ([0, 00) x R'' x T**). We 
may thus apply Lemma 13.31 with cp = if — iip to the last term in (|105p and obtain the following 
asymptotic expansion in ^^([O, 00) x R'^ x T'^), 

e-'^^'''^^''^/'' {{hdt + Op^{a^))e'^/''b){t,x,7j) 

X i{dt^p{t,x,r]) + af{x,iP'^{t,x,ri))) b{t,x,r]-h) + h{{dt + Y + i div y) 5o (?) (i, x, r;) 

00 1^—2 

+ + Y +\dWY)b,^^g){t,x,T^) + Y,{P,^^b.^g){t,x,r^)) , 

where Y and are given by (|95p and p3p with (p — ip — iip. (Here we again use that ip = ip on 
supp (6).) We observe that 

{dt +Y + ldWY){b,g) = {{dt + Y + ^divY)b,) g + m,, 7^ = 0,1,2,... 

PL'->,{b^g) = iP^-><b^) g + m^^t, , < >fr < , = 2, 3, 4, . . . , 

where rn„,fh^^^ G C^([0,oo) x R'' x T"^) are supported in suppdg'j) C {^tp > 0}, which implies 
e»^/''m^ = 0{h°°) and e'^^'^m^^^ = 0(/i°°) in 5b(R x R'^ x T"*) . 

We next use (|89p and the fact that bo,bi, . . . solve the transport equations (To), (^i), ■ • ■ • The 
error terms in (fSQ]) and ([98|) lead to contributions of order 0{h°°) in S'6(R x R'' x T''), too, because 
of the elementary estimates 

□ 
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7 The construction of ^"(0)"^ 

Proposition lG .31 supplies the piece of the paramctrix for E"{Q) corresponding to the set Kq x (27rZ)'' 
as we shall see in the proof of Proposition 17. II below. Outside Kq x (27rZ)'* the Weyl symbol 
is elliptic and we can use standard methods of pseudodifferential calculus to construct the missing 
piece of the parametrix. Here are the details: 

If were elliptic, we had a well-known asymptotic expansion, q{x, ^) x X^l^o 0^ 
the Weyl symbol of the inverse operator (Op)^(a^)) ^ . We can of course always write down this 
expansion formally and determine at all points (x,^) where is invertible. We proceed in 
this way and pick some x € C^{R'^, [0, 1]) such that x = 1 in a small neighbourhood of Kq and 
supp (x) C {x — !}■ We recall that x ha-s been introduced in ([M]) . Then we define the symbol 
q E SbiR"^ X Rrf) to be a Borel resummation of 

oo 

where now each term is a well-defined function. We let denote the composition of Weyl 
symbols, i.e. 

= /'e^<?-''l--)/"+»(e-CI--)/"a(^,,7)c(^,C)^J^^^, (106) 

J yZTrnJ 

for a, c e SbiR"^ X R^); see, e.g., 12, §2.7]. Then we get, by definition of q,y, v e No, 

l-.\a+p\ /I \ \a\ 

a^#r9>,0 ^ E i2^)\'^mL^ (aggfa^(x,0) {dtd^.q{x,0) 

oc 

= l-x{x) + Y,h''r,{x,i), (107) 

i/=0 

where each error term r^, u £ No, contains some partial derivative of x, which shows that supp C 
{x = 1}, for all V e No. We now set q := g#)^(l — x) and infer from p07p that 

in Sbi'R'^ X Rrf). We observe that the symbol q is again 27r-periodic in ^. In fact, if we suppose that a 
and c in p06p are 27r-periodic in the momentum variables, then we see that (a#)^c)(a;, ^4- 27r£) = 
(a#)^c)(a;, ^), £ S Z"*, by a translation of the variables 77 and C, in the oscillatory integral (|106p . 
We let Kq{x,y) denote the distribution kernel of O'p^ [q). The periodicity of q implies that 

K,{x,y) = h^J2 q{^,^)Sz{y-x), x,yeZt, 

where the Fourier coefRcients, g, are defined as in 
Recalling that b vanishes, for large t, we set 

C30 

Iix,y) := J Je^^^''^'^y'^'^^y\^^/Hit,x,v;h)^^^, (108) 

for all X, ?/ e Z^, and try the operator V given by 

Vix.y) := lix.y) + h^q{^,y-x), x, y £ Z^ (109) 
as an approximate inverse for e'''/''£"'(0)e^'^/''. 
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Proposition. 7.1 It holds V e if(£2(Z^)) and 

where TZ e if (^^(ZfJ) such that ||7^|| = OAr(/i^), /or every TV G N. /n particular, 



E"{0)) ^ = e-'^/''P ( ^(-l)'"7^' 



m=0 



Proof: Applying successively Lemma [3.21 Proposition 16.31 and the fact that b vanishes, for 
large t, we get, for all x,y £Zf^, 

(eV/^E"{0)e-^^'' ol){x,y) 

( e^/''-S"(0)e-^/'' e'^it.:v)/h ^(^^ ^ . (^) y | r, >/ft , dt drj 



(27r)^/i 

T'' 

OO 

dt drj 



jd 



( Opr (a^) e»'/'(*.-.'')/'' b{t, ■,rj;h)) {x) e'^^ « I " 



{2nYh 



Jd 



where, for iV G N and x,y G Zf, 



e«(-!.l'<>A^(a;)^ +7^l(a;,y) = x(:c) '5., + 7^l(a;, y) , (110) 



n,{x,y) /^^ / / e^^(''^'^y'^-^(y^^y'^f^{t,x,rj;h)^l^. (Ill) 



T'' 

Since the integrants in (jllip vanish, if t is sufficiently large and fN{t,-,ri;h) is supported in 
+ {\x\ ^ i?}, we have 

max{7^l(x,y) : dist(x,jr) sC i?, y G Z^ } = ©^(/i^) , iV G N . 

Moreover, since 2;, 77) — x\ is uniformly bounded on compact sets, we have |'0^(t, x, rj) — y\ ^ 

\x — y|/2 provided \x — y\ is sufficiently large. Integrating by parts in (|lllll by means of the 
operator 



h{ip'r,{t,x,r]) - y 


V.) 


i 


ij'^{t,x,Tj) - y 


2 



we therefore find that 

|7^l(x,y)| - 0(/i^(a;-y)-^), TV G N . (112) 

If we form an operator, TZi, whose matrix element at (x, y) E Zf^x Zf is given by Tli{x, y), we thus 
obtain an element in .if (^i(Z^))nif(^°°(Z^)). By interpolation we then also have 7^l G ^(^^(Z^)). 
Next, we notice that 

Oplfa-X + O -Opr(a^)oOpr(g) = Op,(a^#,g), 

where r := a^#f q - (1 - x) is of order 0{h°^) in S'6(IR'' x R^), 

dyd^ 



{2TThY ' 
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and Op^(a^#/i(7) denotes the /i-pseudodifferential operator defined by the symbol Oy^.jj^hq G 



Sb{R'^'^ X Rjj). The symbols r and a'^^hq arc both 27r-periodic in the momentum variable. More- 



for x,y & Zf, where the Fourier coefficients, {d^^hq)^{x,y, z), z e Z^, are given by a formula 
similar to (|21|) . Since any pseudodifferential operator determines uniquely its distribution kernel, 
we obtain 

(e^/''S"(0)e-^/'' h'^ q{^,y - •)) {x) 



= {l-x{x))5,y + h''9{^,y~x). (113) 

Setting 

(7^2/)(x) := xezf,,/er(z^), (114) 

where p e {l,oo}, we again have 7^2 G ^(^^(ZfJ) n Jf(£°°(Z^)), because 

which is easily verified. By interpolation we see that || 7?.2 \\se{e'(j_dj'^ — 0{h°°). In view of piop - 
(|114p the proof is complete with TZ :=TZi + Tl2- Q 



8 Calculation of the leading asymptotics 

In this section we calculate the leading asymptotics of 'P{x*,y*), where x* and y* are our two 
distinguished points which satisfy Hypothesis 14.41 Integrating by parts repeatedly in the integral 
defining the Fourier coefficient g( ^ 1"^ , y* — 



r{x\y*) 



X*) we see that 
I{x\y*) + 0{h^) 



(115) 



and it remains to determine the leading asymtotics of I{x*,y*). To this end we recall that H — 
T — U and that Xh denotes the Hamiltonian vector field of if. tt : T*IR'' R'' is the canonical 
projection in the following. We let p : IR ^ r*IR'^ denote the unique integral curve of Xh running 
in f = {H — 0} such that 7 := t^p\[q^t] is a minimizing geodesic from 7(0) = y* to 7(r) = a;*, for 
some suitable r > 0. Recalling the notation introduced in the paragraph preceeding Theorem l2.6l 
we have 

vy. = 7(0) = -k{T,x\Q), (116) 
V,. = 7(t) = Q(t,2;*,0). (117) 

Furthermore, we let X(i), P(t) G .^(R"*), t G R, denote the solution of the initial value problem 



d /X{t) 
dt \P(t) 



X{t) 

Pit) 



XiO) 
P(0) 



(118) 



Proposition. 8.1 The following formula holds, as h\, 0; 



Iix\yn = 



Jd 



dr] dt 
{2t:Y h 



(- 

V27r 



n \ -5- 



X{t) 



+ 0{h—) 



(119) 
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Proof: We apply the method of complex stationary phase [14, §2] (with respect to t and 77) 
to the integral pi9p . The critical points of the phase are given by 



= dttp{t,x\fj), 

= \/^^/j{t,x\7j) - y* . 



(120) 
(121) 



Of course, a critical point gives a non-vanishing contribution to the asymptotics of (jll9p only if 
it is contained in the set = 0}. By CoroUarv 15.101 we know that we must have t = or 
{t,x*,ri) e ^ in this case. At t = we have, however, V,,V'(0, a;*, 77) = x* ^ y* . So, (^0]) imphes 
that (t, 0) is the only critical point of the phase in (|119p where ^'0 vanishes. Again by 



a,V(T,x*,o) 



0. 



To study the remaining second derivatives of the phase we differentiate the identity Q{t^ k{t, x* , 77), rj) — 
X* with resect to t and 77 and obtain 



Q'y{t,k{t,x*,rj),r])k{t,x*,7]) + Q'y{t,k{t,x* ,T]),r]) k{t,x* ,1^) 
Q'y{t,k{t,x*,T]),r])k'^{t,x*,i]) + Q'y{t,k{t,x*,'q),T]) k'^{t,x*,T]) 

Evaluated at {t, rj) — (r, 0) this yields 



■'Q{t,k{t,x*,r]),ri) , 

^Q'rj{t,k{t, X*,l]),f]) . 



Q'y{T,y\0) Vy. 

YyiT, y*,0) i7/^^'^(T,X,0) 



Vx* , 

■iQL{T, y*,0) 



(122) 
(123) 



if we take (fM|) . (|116p . and (|117p into account. We consider iQ',^{T,y* ,0) more closely in the 

following. To this end, we recall that the derivative of the flow of Xa-p along {j{t),0), kJ(j/*,0), 
t Cz [0, r], is a solution of the Jacobi equation 



(124) 



where Fa(7(i),0) is given by (|69p with x{t) = j{t) and e = 0. Since we want to determine 
iQ'rjir, y*, 0) = i 7r^(Kr)J,(2/*, 0), we consider the solution of ()124p with initial condition i(Ko)J,(j/*, 0) 

(0, it), which we denote by {X{t),P{t)). It follows that = and nP = and that (nx, 3P) 
is a solution of the Jacobi equation 



d_mx{t)\ 



-*'E{t) 



(125) 



The matrix in l|125p is the fundamental matrix at {x,p) = (jit), 0) of the Hamiltonian Hip{x,p) := 
H{x,p + if'{x)) = {H oE){x,p), where E{x,p) := {x , p + (p' {x)) , {x,p) G T*IR'*, defines a symplecto- 
morphism. We further recall that H'^ ,^^{x,0) = 0, for a; e Xq, by construction of (p. Using these 
remarks we check that ^X = X, where (X, P) is a solution of (jllSp . In fact, the Jacobi equation 
(fna is solved by {E')-^{-/){X, P) = (X, -(^"(7)^ + P). 

In summary, we see that iQ'^{,T,y* ,0) — X(t). In particular, iQ'^{,T,y* ,0) is real and so is 
iV'(,'^(T, x, 0) by p23p . because Q'y{T,y* ,0) is a real invertible matrix. 

By Lemma [Owe know that ba{T,x*,0)-^ = |Q^(t, y*, O)!^/^. Using also (fH^ and (fHSj) we 
thus get 



bo{T,x*,0) 







2 


1 


2 




k'.iy*) 








-*7(0) 


2 


-'Vy^ 


7(0) 






v,. X{t) 



where all derivatives of -0 are evaluated at (r, a;*',0). From Proposition 18.21 below it follows in 
particular, that the last determinant in the previous equation is non-zero. Therefore, the left 
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hand side is non-zero, too. (To show that i/j'pj (r, a;* , 0) is non-singular we could also appeal to 
the strict positivity of Ct at (x*, 0, y*, 0).) We may thus split the integral pi9|) into pieces by 
means of suitable cut-off functions and apply the stationary (complex) phase formula T31 §2] to a 
sufficiently small piece of (|119p near the critical point (r, 0). 

Finally, we remark that, for all (5,T > 0, inf{!3V(i, a;*, 77) : t e [5,T], rj e J'^, ^ S} > and 
the corresponding piece of the integral (|119p gives a contribution of order 0{h°"). At i = we 
have V^(V'(0, x* , ri) — {y* \ r])) = x* — y* ^ 0, rj £ J"^, and, for 77 = and some S' , 5" > 0, 

\V^\^=omt, x*,7])-{y*\v))\^ m, x\0)-y*\^5', t e [5,t - S" ]U [t + d",T] , 

Therefore, the remaining pieces of the integral (jll9p can be treated by means of integration by 
parts in rj and give contributions of order 0{h°°), too. □ 

It is of course simple to check that (|119p agrees with the already known formula ITCl) in the 
translation invariant case. In the next proposition we obtain a generalization of (|17p . 

To prove it we introduce special orthonormal bases with respect to the Finsler structure O 
pp. 31] in the following. 

We denote the quadratic form associated with G{z,Vz) by gz, for z — x*,y*, and pick a 
basis, (61,..., which is orthonormal with respect to gy* such that bd — F{y* ,Vy*)~^Vy* . 
Similarly, we pick another basis, ( ci, . . . , ), which is orthonormal with respect to gx* and where 
Cd — F{x* ,Vx*)~^Vx* ■ We denote the dual basis vectors by b*, c*, 1 ^ i ^ d, so that b*{bj) — 
= etc. 

For l^i^d— l,wc further introduce a Jacobi field, Ji, which is by definition the solution of 
the Jacobi equation 

D^DsJr - R{J^,^)^, MO) = b,. 

Here ^ is the velocity vector field of the unit speed geodesic from y* to x* , denotes covariant 
differentiation in the direction ^ with reference vector and 

d 

k^£,m,n—l 

where Rk'mn are the components of the /i/i-curvature tensor of the Chern connection defined by 
F [5] . Since we only need a well-known Taylor expansion of mutual scalar products of the Ji and 
the formula 

(eyiY>y)'{rbd)b^ = - J,(r), 

where 

r := dF{x*,y*), 

we do not explain these notions further. 



Proposition. 8.2 The following identity holds: 





Vx* X{t) 



G{x\vx*)G{y\vy.) 



1 1/4 



\/{Px* \ Vx* ){Pr \ Vy 



[gx*iJiir),Jj{r)) 



1/4 



whe 



{9x*{Mr),Mr)))'' ' ^ dF{x\y*)^ + 0{dAx\y*)^) 



1/4 



as X — > y 
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Proof: To simplify our notation we denote the distinguished points x* and y* simply by x 
and y in this proof. We let B and C denote the matrices whose i-th row is h* and c* , respectively. 



From = F'^{x,Vx) = Px we infer that C Vx 



'iO...OF{x,Vx)) - {px\vxyi0...01) and 



similarly Vy B 



1 
C 



|uy)(0...01). We thus get 



— Uy 

X{t) 



1 

*B 



[Px\Vx){Py\Vy] 



CX{t)'B 



Here the determinant on the right side clearly equals 



det 



The vectors 1 ^ j ^ d — 1, are orthogonal to bd G Sj, with respect to the Euclidean scalar 
product and, hence, span the tangent space of the figuratrix at the momentum which is conjugate 
to bd-, i.e. at py. It therefore suffices to consider 

^iT)\Tpy<Sy^ i^T\Sy)'{Py) , where a;^(p) := 7rexp(rXff)(y,p) . 

Let C : ?z ^ denote Legendre transformation, so that C maps p g (Fz to its conjugate direction 
V e 0^, for z e R''. It holds C{px) ~ Cd \\ Vx and C{py) — bd \\ Vy. Moreover, we have p. 410] 



C'iPy) = G{y,C{Pv))-'\u^ G{y,Vy)-'\ 



T (F 



(126) 



We recall that G is homogenous of degree zero in v. Besides we know that 
((exp )'(r£(p.y)) £(p.y) = exp {tC{py)) || Vx , 

az t—r 

where r — dp{x,y). The image of C{py) under (expj^)'( r£(py)) is therefore -orthogonal to Ci, 
1 ^ « ^ d — 1. (This is precisely the statement of the GauB lemma in Finsler geometry §6.1. A], 
as the vectors c,, l^i^d— 1, span ^^(p^jSi:.) In particular, 

c* (exp^)'(r/:(p,)) [C{py) ® (dpixri-), y))'] - . (127) 

We shall use the following identity, 

Xr{p)=expy(dFixr{p),y)C{p)), pe^y. (128) 

Differentiating (fT^ and using P?7| . and c*(u^) = 0, 1 s$ i s$ d - 1, we obtain 

cn^rhJ{Py)'b* = rc*{iexpJ{rC{py))Giy,Vyyb* 
= rc* {{expyy{rbd))bj , 

for 1 ^ i, j ^ d — 1. Let M be the matrix with entries c* ((expj^)'(r 6^)) bj, 1 ^ i, j ^ d — 1. It 
holds J2i=i *c* c* = G{x, Ua;)['T£(p^)S,, whcucc we get 

(*MM)^^. = .94(exp^)'(r6,)6,;, (exp^)'(r6,)6,) 

2 

= - y 9v{RiK bd)bd , b,) + Oir^) . 

Here the well-known Taylor expansion in the last line is derived using the formulas of [SJ §5.5]. 
We conclude the proof of the proposition by noticing that (dcti?)^ = dct(*i?i3) = det G{y,Vy) 
anddetC= (detG(x,w^))i/2. ' □ 
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A Applicability of a result by Bach and M0ller 

In this appendix we show that we can apply the resuhs of ^ in our situation and thus obtain a 
proof of Theorem 12.41 The starting point will be a Helffer-Sjostrand formula established in [13] . 
It involves the inverse of a certain operator which is introduced in the following. First, we recall 
that we have fixed some pure, tempered Gibbs measure, fJ,f3^h, and have set Jf'^ = L'^{^j.p,h)- We 
further set Jff^ := L^(^/3^/i)(8)i'^(Z^), where (g) deontes the (completed) tensor product of Hilbert 
spaces, and 

{ / : R^' ^ C : 3 n e N , XI, . . . ,x„ e /„ e C^{R") 
Vae R^' : /(a) = /„K,, . . . , } . 

Then 17" is a dense subspace of and the algebraic tensor product H.^ := (g) {cx ■ x € ) is 
a dense subspace of . Here {cx ■ x G Zf) is the linear hull of the canonical orthonormal basis 
of i'^i'I.f)- We next define two operators, -S^'^*),, « = 1, 2, on with domain dom(^^*]^) :— by 

for cr e R^^ and / G 17°, and 

We recall that the matrix element of E" at ct £ R^'' is defined by (fTTj) . Then -Sf^*)j is essentially 

selfadjoint on J^*, for i = 0, 1 [131. We denote its unique selfadjoint extension by We recall 

a result from [13] which permits to locate spectral gaps at the bottom of the spectrum of ,j at 
large inverse temperatures. To this end, we set 

:= inf i?^'e(x,0) - J V sup l^^'e(a;, £, 0) , 

:= sup D'/,g{x, 0) + 2 J ^ sup I, 0) , 

and assume that there is some n g N such that 

(n + l)m^ - TITO* > 0. (129) 

Note that this is always true, for fixed n, provided J > is sufficiently small. Moreover, we set, 
for £ > 0, 

n 

J{n,e) := \^ {v rrii, — e , v m* + e) . 

Theorem. A.l For every e > 0, there is some /3o 1 such that, for all (3 ^ /3o and h G (0, 1], 
spec(i^°)^) n (-CX), (n+ 1)to* -e] C I{n,e)U{0}, 
spec(L^^^)j) n (-00 , (n + 1) TO* - e] C T{n,e). 



Next, we introduce an effective Hamiltonian, or, Feshbach operator, Tp^h, associated with i^^jj. 
To this end we denote the projection on onto the subspace of constant functions by p, and 
set P := p i8) 1. We write p^ = 1 — p and = p-*- (g) 1. In view of Theorem I A . 1 1 !Fr h is then 
well-defined by 

:F,.. ■= PL^iP-PL^iP^{P^L^iPT'P^LfiP 

= PE" P - PE" P^{P^ L^p^P^y^P^ E" P . (130) 
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Now we are in a position to write down the Helffer-Sjostrand formula for our lattice spin model. 
Its validity is also proved in [13]. We note that it only holds if we assume to be pure. The 
crucial point is that kerLjj'l^ C 1 if and only if ^pji is pure. 

Theorem. A. 2 There is some Po ^ 1 such that, for all (3 ^ /3o, h G (0, 1], and x,y £ if^, 

Cot p,h{(J^\ (jy) ^'^{e^ \ {L^^^)'^ Cy) = i^P^hY^^y) ^ (131) 

where ez = 1 ^ Cz, for z E T'l- 



To obtain a proof of Theorem 12.41 we only have to verify certain conditions that allow to apply 
[H Theorems 3.3 & 3.7] to our situation. In fact, it suffices to derive a certain estimate involving 

(i) . . . 

the operators /j which is done in the next lemma. The main idea behind all this is that, due 
to the small temperature localization at the global minimum of the spin system, we expect the 
Feshbach operator J^/s^h to be well-approximated by E"{0); see (|130p . Before we come to that we 
have, however, to introduce some further notation. 

We let De e ^{P{Z'l)) denote the diagonal part of E"{0), that is, Dec^ := E"{0)^^e:c, for 
x G Z^. Furthermore, we set 

T := D-^'^^{De -J-'E" iO))D-^'/\ 
V := (1®D£'/')(S"-10^"(O))(1®D^'/'), 
Bo := ijj"' ® + 1 , 

B lJj"^ (8)0^1 + 1® (D^^/^£;"(0)D£^/^) = Bo - J1«)T. 

Identifying ^^(Z^) with Ran P = CI ® ^^(Z^) we have 

(1 ® D^'/') J^f3,h (1 ® D^'/') = 1 - JT + Y , 

where 

Y := /3i/2pvp - p^/^pyp^{p^{B + y)p^y^p^yp. 

Lemma. A. 3 There exist ^ 1, Jo > 0, and d E (0, 1) such that, for all (3 ^ /3o, J G (0, Jq], 
and x,y E J-f^, 

\y.v\ ^ {(i-^?jT)-i}^^. 

Proof: Combining the estimate succeeding Equation (53) of [13] and Lemma 7.4 of [13] we 
see that there is some constant C E (0, 00) such that, for all sufficiently small J > and all 

{E'iMKyi^)y"^ \Ky - KyiO)\ ^ ^ i5.y + JT.y) {Lf + 1) (132) 

in the sense of quadratic forms on Vf' . Here we also use that we have a strictly positive uniform 
lower bound for i?"^(0), x E T'l, provided J > is sufficiently small. Moreover, we use that there 
is some c > such that, for all x,y E Z^, we have the bound J^j, ^ cJT^y, where the Jxy ^ 
are the constants appearing in [TJl §7]. The form bound (|132p clearly implies that 

\{e \PyPe ) I - KM(-g"y--g"y(0))l),^o| < _^(s ,jj ) 

^E'^ME'^yW ^' 

because L^^^ 1 = 0. Since 

oo 
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due to the fact that T^y ^ with strict inequahty at least for \x — y\/h = 1, there is some i? € (0, 1) 
such that 

/3i/2|(e.|PVPe,)^,| C{(l-^JT)-i}^^. (133) 

We set Bo := P-^BoP-^, B := P^BP-^, and V := P^MP^ in the foUowing. Since the numbers 
E'^^iO) are bounded from above uniformly in x £ Zf , it also follows from (|132p that there is some 
C £ (0, oo) such that 

II {B-'^^yB-'/%y II + II {Bo '/'V Bo II ^ ^ {S,y + JT^y) (134) 

for all x,y £ Z^. Now let e £ (0, ni* — m*). By Theorem l2.4f i') we know that, for sufficiently large 

— rUi, — £ 1 
Bo ^ — ^ + 1 =: 



where £ This implies, for all x, y e ZJ^, 

II {Bq ^ "^TBo ^ la;^ II ^ >^T2;y||BQ II ^ d' JJxy, 

which permits to get 

oo 

||{Bi/'B-'By'}.,|| ^ ||{^(Bo^/'jTBo'/y} ^ {{1- Jjy^U . (135) 

M ^ ^ — ' J xy 

n— 

Using (|134|) and (|135p we can now copy the proof of [H Theorem 4.6] to get 

/3i/2|^e,|PVP^(P^(B + V)P^)"'p^VPe,y)| C {{1 - d' Jl)~^}xy , 

for some C" € (0, oo). The previous estimate together with (|133p implies the assertion of the 
lemma. □ 

The proof of Theorem l2.4f iii') now follows from Lemma fA. 31 together with f?, Theorems 3.3 & 3.7]. 

B Some elementary linear algebra 

In this appendix we prove that Formula (fT6| of Theorem 12.51 yields Formula (fT7|) . 

We let iV : (F ^ S'^~^ denote the exterior normal field on IF. (Note that, of course, the 
figuratrix does not depend on x in the translation invariant case.) The derivative of its inverse 
is given by {N'^^y{v) = \Hp{p{v))\ (_ffpp(p(w))^) , for v £ S'^^^ . On the other hand we have 
F{v) = {v\N~^{v)) and direct calculations show that the Hessian of P at u restricted to the 
orthogonal complement of v in IR is given by P^;(w)-L = ^ (7V-i)'(^). (Here and in the following 

we always assume w e R"^ \ {0} and set v :— v/\v\.) This implies 



d-l 



*?; 



(136) 



Next, we pick a basis, {bi, . . . ,bd), which is orthonormal with respect to G{v) — \ (P^)";^^ (w) and 
such that bd — F{v)^^v. We denote its dual basis by (6*, . . . jb'^). Furthermore, we denote the 
matrix whose i-th row is b* by B and the matrix whose i-th column is 6; by Q, so that B Q — Id- 
Then [51 J 2. 2] 

F:,iv) = 'B( T^^'-' Mp, \Q\-' = \G{v)\, CQv)d= 



; ' ' ' ' F{v) 
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Consequently, we get 



*w 



\Q\- 











^G{v)\ 



(137) 



From (|136p . (|137p and the liomogenity of F we infer that 



d-3 



\G{v)\ 



mp{i))^\\v\''~' \H'^(p{vWF{vY F{v 



Finally, using the fact that N ^{v) — p{v) and v \\ WpH{p{v)) we see that 

\H;ipiv))fFivr = {V,Hipiv))\piv))\ 
This shows that (fTH]) and (fT7|) are equivalent. 

C Basic facts about almost analytic extensions 

We recall some notions developed in [H] : always denotes some open subset of C" in the following 
and riiR := f2 n R". Two functions, 31,52 & C°°{n,C), are called equivalent, in symbols gi ^ 32, 
iff, for every compact il' (s and G N, there exists a constant CM,n' G (Oj °°) such that 

We recall a technical lemma from |14) which is used very often in the main text. 

Lemma. C.l Let Q C R" be open, f e C^{Q), and g : Q ^ R be locally Lipschitz continuous, 
i.e. for all compact Q' ({= Q there is some Cg/ € (0,oo) such that \g{x) — g{y)\ ^ Cgi\x — y\, 
x,y (z Q' . Assume that supp [g) = Q and that, for all compact Q' is= Q and N G N . there is some 
Cn g' G (0,00) such that 

\f{^)\ ^ CN,g'\gix)\'' , xeg'. 

Then, for all compact Q' <s= Q , N Cz N, and a G Nq, |a| < N, there exists some C'N,g',a G (0,oo) 
such that 

\dSf{x)\ ^ CN,g',a\gix)f~\''K xeg'. 



Let gi,g2 G C°°{il, C) be equivalent. By the previous lemma we find, for every compact il' d il, 
G N, and all multi-indices a,/3 G Nq, some CN,n',a.i3 G (0,oo) such that 



z G SI'. 



A function / G C°°(ri, C) is called almost analytic, iff dg-f ^ 0, for i = 1, . . . ,n. We denote the 
set of all almost analytic functions on fl by C'*(17,C). Furthermore, let / G C°°(17|r,C). Then 
/ G C°°{Q,C) is called an almost analytic extension of /, iff / G C"^''(ri, C) and f\n^= f. We 
denote the set of all almost analytic extensions of / by C^^[f]. For instance, every / G C^(IR", C) 
possesses almost analytic extensions and /i, /2 G C'^'^[f] fulfill /i ^ /2. 

Conversely, if / G C^'^[/] and g ^ 0, then again f + g & C'^^lf], so that C^^[f] is an equivalence 
class modulo ~. For every S > 0, there is some / G C"^'^[/] with support contained in {z G 
il : 5}. Moreover, for / G C'^if], for aU multi-indices a,/3 G Ng, and i = l,...,n, 

it holds d^z^Qz'^if ~ by Lemma [C.ll In particular, it follows that any partial derivative of 
/ G C^^{n, C) is again almost analytic. 
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